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SUMMARY 


This  report  is  the  technical  report  on  the  second  phase  of 
an  investigation  of  the  hinge  moments  of  all-movable  controls  as 
utilized  on  cruciform  missiles.  In  the  first  phase,  an  attempt 
was  made  to  develop  a  hinge-moment  prediction  method  as  general 
as  possible  for  speeds  ranging  from  subsonic  to  hypersonic.  In 
that  work,  the  method  was  successfully  applied  for  Mach  numbers 
ranging  from  1.3  to  3.7.  In  the  present  phase  of  the  investiga¬ 
tion,  the  methods  from  phase  one,  with  some  modifications,  have 
been  used  to  develop  a  procedure  for  determining  the  optimum  fin 
planform  and  airfoil  section  to  minimize  the  fin  hinge  moment 
over  a  range  of  Mach  numbers,  angles  of  attack,  and  fin  deflec¬ 
tion  angles.  The  Mach  number  range  is  restricted  to  supersonic. 
The  transonic  speed  range  cannot  be  handled  since  methods  are 
not  presently  available  for  predicting  the  center-of-pressure 
shift  due  to  thickness  in  this  speed  range.  Experimental  data 
are  not  available  to  assess  the  accuracy  of  the  method  at  hyper¬ 
sonic  speeds. 
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PREFACE 


This  technical  report  covers  the  work  performed  under 
Contract  N00014-81-C-0276  from  December  1,  1981  to  December  31, 
1982.  The  program  is  sponsored  by  the  Office  of  Naval  Research, 
Arlington,  Virginia.  Dr.  Robert  E.  Whitehead  is  the  Scientific 
Officer. 

Support  for  the  program  was  provided  by  Air  Force  Wright 
Aeronautical  Laboratory  (AFWAL) ,  under  the  leadership  of 
Mr.  Vernon  0.  Hoehne. 
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1 .  INTRODUCTION 


The  determination  of  the  fin  size,  planform,  and  airfoil 
section  to  provide  a  missile  with  adequate  control  while  minimiz¬ 
ing  fin  hinge  moment  has  received  little  analytical  treatment 
over  the  years.  The  problem  has  normally  been  studied  by  per¬ 
forming  expensive  and  time  consuming  wind-tunnel  testing.  Since 
the  weights  of  the  control  actuators  and  their  power  supplies  are 
substantial,  the  design  of  a  minimum  hinge-moment  fin  would  help 
keep  this  weight  down. 

The  work  described  in  Reference  1  was  directed  toward 
developing  a  method  which  would  accurately  predict,  including 
thickness  effects,  the  hinge  moments  on  all-movable  controls. 

It  uses  the  methodology  developed  in  References  2  and  3  and  ex¬ 
panded  upon  and  modified  the  computer  program  of  Reference  3, 
MISSILE2.  It  was  found,  in  Reference  1,  that  the  canard  hinge 
moments  were  well  predicted  at  supersonic  speeds  except  for  lee¬ 
ward  fins  when  strong  body  vortices  were  present.  At  transonic 
speeds,  the  predictions  were  not  good.  As  a  result  of  these 
comparisons,  it  is  felt  that  the  methodology  exists  for  develop¬ 
ing  a  method  for  determining  the  optimum  fin  planform  and  airfoil 
section  for  minimizing  the  hinge  moment  acting  on  all-movable 
controls  at  supersonic  speeds.  It  is  the  purpose  of  this  report 
to  summarize  such  a  method. 

This  work  is  being  conducted  under  contract  N00014-81-C-0267 
from  the  Office  of  Naval  Research  which  is  supported  by  the  Flight 
Dynamics  Laboratory  of  the  Air  Force  Wright  Aeronautical  Labora¬ 
tories.  This  is  Phase  II  of  the  work.  The  results  of  Phase  I 
are  documented  in  Reference  1. 

In  this  report  a  general  discussion  of  the  problem  is  first 
presented.  This  is  followed  by  an  explanation  of  the  modifica¬ 
tions  made  to  program  MISSILE2  and  an  example  of  how  the  optimum 
hinge-line  location  is  determined  for  a  particular  fin.  Next,  the 
procedure  used  to  determine  the  optimum  fin  to  minimize  hinge 
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moment  for  a  particular  set  of  constraints  is  described.  The 
last  section  presents  some  concluding  remarks.  The  computer  pro¬ 
gram  is  described  in  an  appendix. 

2.  GENERAL  CONSIDERATIONS 
2.1  Reasons  for  Study 

The  use  of  all-movable  controls  on  missiles  for  trim  and 
maneuvering  is  quite  common  since  the  greatest  effectiveness  is 
very  often  obtained  by  using  all  of  the  fin  area.  Since  missiles 
operate  over  a  wide  Mach  number  range,  the  longitudinal  center- 
of-pressure  location  of  the  force  produced  by  the  fin  can  vary 
substantially.  This  variation  can  be  increased,  further,  by  the 
requirement  of  high  angle  of  attack  operation  for  greater  maneuver¬ 
ability  and,  thus,  the  requirement  for  larger  control  deflection 
angles.  At  the  higher  angles  of  attack  and  control  deflection 
angles,  nonlinearities  are  introduced  which  make  the  prediction 
of  the  center-of-pressure  location  very  difficult.  The  accurate 
prediction  of  this  location  over  the  range  of  flight  conditions  is 
important  since  this  location  relative  to  the  control  hinge-line 
position  affects  the  hinge  moment  acting  on  the  fin  and  this,  in 
turn,  affects  the  capabilities  required  of  the  control  actuators. 
Greater  capabilities  increase  the  weight  of  the  actuators  and 
their  power  supplies.  In  order  to  minimize  this  weight,  it  is 
important  that  an  optimum  fin  planform  and  airfoil  section  be 
selected,  in  addition  to  the  control  hinge-line  location,  so  that 
the  hinge  moment  is  kept  to  a  minimum  over  the  range  of  flight 
conditions  expected. 

Until  recently  preliminary  design  methods  have  been  available 
only  for  the  linear  range  of  angle  of  attack  and  control  deflec¬ 
tion  for  predicting  control  normal  force  and  hinge  moments.  In 
Reference  1  a  method  was  presented  for  estimating  control  hinge 
moments  to  large  angles  of  attack  and  control  deflections.  It 
was  found  there  that  the  predictions,  at  supersonic  speeds,  were 
quite  accurate  as  long  as  strong  nose  vortex  effects  were  not 
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present  and  the  forward  fins  were  not  producing  strong  interfer¬ 
ence  effects  on  the  rearward  fins.  At  transonic  speeds,  all 
predictions  were  inaccurate.  As  a  result  of  this  work,  the 
present  investigation  was  undertaken  with  the  goal  of  developing 
a  method  for  determining  the  optimum  fin  planform,  airfoil  sec¬ 
tion,  and  hinge-line  location  for  a  single  set  of  cruciform  fins 
on  a  body  of  revolution  at  supersonic  speeds. 

2.2  Optimization  Factors 

A  number  of  factors  may  affect  the  determination  of  the 
optimum  fin  planform  and  airfoil  section.  One  of  these  is  the 
range  of  flight  conditions  which  the  missile  will  encounter  dur¬ 
ing  its  trajectory.  These  conditions  are: 

Mach  number,  M 

'  oo 

Dynamic  pressure, 

Combined  angle  of  attack,  ac 

Roll  angle,  cf> 

Fin  deflection  angle,  6 

The  complete  range  of  these  parameters  should  be  investigated 
since  they  each  individually  and  in  combination  affect  the  fin 
hinge  moment. 

Another  factor  which  may  influence  the  selection  of  the 
optimum  fin  is  the  force  which  the  fin  produces.  This  force  must 
be  sufficient  to  give  the  missile  adequate  trim  and  maneuver¬ 
ability.  Thus,  certain  fins  may  be  ruled  out  because  they  do  not 
produce  sufficient  control. 

A  factor  which  could  reduce  the  number  of  fin  planforms 
which  can  be  considered  would  be  a  restriction  on  fin  span.  Be¬ 
cause  of  the  way  the  missile  is  carried  or  launched,  a  maximum 
span  restriction  may  be  imposed.  There  may  also  be  restrictions 
applied  to  fin  taper  ratio  or  thickness  for  structural  reasons. 

The  above  are  some  of  the  factors  which  may  restrict  the 
ranges  of  some  of  the  variables  which  must  be  considered  in  an 
optimization  study.  In  the  present  work  and  computer  program  the 
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fins  are  restricted  to  unswept  trailing  edges  and  the  planform  is 
characterized  by 

Fin  exposed  semispan,  (s  -  a) 

Fin  root  chord,  c 
'  r 

Fin  taper  ratio,  A  =  c,/c 

t  r 

From  these,  the  aspect  ratio  and  fin  area  are  determined. 

Fin  aspect  ratio,  JR  =  ^  ^ 

r 

c 

Fin  area,  S_  =  (1  +  A) (s  -  a) 

The  fin  aspect  ratio  is  defined  to  be  that  of  the  fin  and  its 
mirror  image  joined  together  at  their  root  chords.  Three  of  the 
above  five  parameters  must  be  specified  and  varied  in  the  opti¬ 
mization  study.  The  fin  airfoil  section  assumed  in  the  computer 
program  is  shown  in  the  following  sketch. 


In  sections  parallel  to  the  root  chord,  it  is  assumed  to  be 
similar  across  the  span  of  the  fin  and  consists  of  a  flat  plate 
of  length  f  with  equal  length  wedges  for  the  leading  and  trailing 
edges.  For  the  purpose  of  an  optimization  study,  it  is  charac¬ 
terized  by  two  parameters. 

Thickness  to  chord  ratio,  t/c 
Flat  section  to  chord  ratio,  f/c 
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There  are,  thus,  five  fin  parameters,  three  for  the  planform  and 
two  for  the  airfoil  section,  which  must  be  varied. 

2.3  Use  of  Program  MISSILE2 

The  computer  program  MISSILE2,  Reference  3,  is  an  engineer¬ 
ing  prediction  method  for  determining  the  forces  and  moments  on 
cruciform  missiles  to  high  angles  of  attack.  It  predicts  fin 
normal  force  quite  accurately  but  does  not  predict  the  fin  axial 
center-of-pressure  location  with  sufficient  accuracy  to  give  a 
good  hinge-moment  prediction.  In  the  work  of  Reference  1, 

MISSILE2  was  modified  in  order  to  improve  the  center-of-pressure 
prediction  and,  hence,  the  hinge  moment.  The  modifications  made 
were 

1.  Provide  accurate  wing-alone  normal  force  and  center-of- 
pressure  positions  as  a  function  of  angle  of  attack  as 
input. 

2.  Add  the  capability  for  determining  the  effect  of  free 
vortices  on  fin  axial  center-of-pressure  position. 

3.  Add  a  better  means  for  extrapolation  outside  the  A* =  2 
limit  of  the  data  base  of  MISSILE2. 

4 .  Add  a  better  means  of  extrapolation  above  Mm  =  3  in  the 
data  base. 

5.  Change  the  method  of  accounting  for  fin-fin  interference. 

All  of  these  modifications  are  described  in  section  4  of 
Reference  1.  With  the  exception  of  modification  1,  they  have  all 
been  retained  in  the  present  program.  Under  modification  1,  the 
wing-alone  normal-force  coefficient  and  center-of-pressure  location 
were  tabular  input  to  the  computer  program.  The  present  version 
uses  the  wing-alone  data  base  contained  in  MISSILE2  to  determine 
the  wing-alone  normal-force  coefficient.  MISSILE2  has  been  modi¬ 
fied  to  calculate  the  wing-alone  center-of-pressure  location 
including  thickness  effects  as  a  function  of  angle  of  attack  there¬ 
fore  eliminating  that  as  input  data.  The  methods  used  are  described 
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in  section  4.2  of  Reference  1.  Some  other  modifications  have 
also  been  made.  All  modifications  will  be  described  in  the  next 
section . 


3.  MODIFICATIONS  TO  PROGRAM  MISSILE2 

A  number  of  modifications  have  been  made  to  program  MISSILE2 
since  the  work  described  in  Reference  1.  They  are 

1.  Specification  of  the  fin  airfoil  section  as  input  data 

2.  Calculation  of  the  flat-plate  axial  center-of-pressure 
location  for  the  fin  (no  thickness  effects) 

3.  Calculation  of  the  axial  center-of-pressure  shift  due 
to  fin  thickness  distribution. 

4.  Determination  of  fin  axial  center-of-pressure  location 
at  high  angles  of  attack 

5.  Tabulation  of  maximum  and  minimum  hinge  moments  as  a 
function  of  Mach  number  and  hinge-line  location 

In  the  following  sections,  the  above  items  and  their  implementa¬ 
tion  in  the  computer  program  will  be  described.  A  full  descrip¬ 
tion  of  MISSILE2  will  not  be  given  since  this  information  is 
contained  in  References  2  and  3.  Similarly,  the  modifications 
made  during  the  work  of  Reference  1  will  not  be  described  since 
that  reference  contains  a  complete  description. 

3.1  Airfoil  Section  Specification 

The  fin  airfoil  section  is  used  in  determining  the  axial 
center-of-pressure  shift  due  to  fin  thickness.  This  calculation 
requires  knowledge  of  the  local  surface  slope  of  the  airfoil  as 
a  function  of  distance  along  the  chord  at  a  series  of  locations 
across  the  span  of  the  fin.  In  the  modification  presently  in 
the  computer  program  the  assumptions  have  been  made  that  the 
airfoil  section  is  a  flat  plate  with  equal  wedge  angles  in  the 
streamwise  direction  for  the  leading  and  trailing  edges  and 
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that  the  airfoil  sections  are  similar  across  the  span  of  the  fin. 
The  airfoil  section  is  shown  in  Figure  1  and  is  defined  by  two 
parameters;  the  thickness  to  chord  ratio,  t/c ,  and  the  flat-plate 
length  to  chord  ratio,  f/c.  These  two  quantities  are  input  to 
the  computer  program  as  is  the  number  of  strips  across  the  span 
of  the  fin  to  be  used  in  the  strip-theory  calculation  of  the 
center-of-pressure  shift  due  to  thickness. 


3.2  Axial  Center-of-Pressure  Location 
Without  Thickness  Correction 

The  axial  center-of-pressure  location  for  a  fin  with  no 
thickness  is  determined  using  linear  theory.  This  quantity,  for 
most  planforms,  can  be  obtained  from  figure  4.1.4.2-26  of  DATCOM, 
Reference  4.  For  delta  wings  the  center  of  pressure  is  always 
at  the  2/3  root-chord  location.  For  rectangular  wings  with  effec¬ 
tive  aspect  ratios,  37R,  equal  to  or  greater  than  one 


nt 


lfl  -  2/ (3  ml 

2  |_1  -  l/(2$m 


(1) 


where  3 


= 


M  -  1  and  the  aspect  ratio,  7R,  is  that  of  two  fins 


joined  together  at  their  root  chords, 
x/c^  is  obtained  from  tl 
Chart  10  of  Reference  5, 


If  3®  is  less  than  one. 


x/cr  is  obtained  from  the  following  table  which  was  obtained  from 


3A* 

0 

0.5 

1.0 


x/c 


nt 


0 

0.2 

0.333 


The  determination  of  the  center  of  pressure  for  the  fin  with 

no  thickness  using  the  above  sources  of  information  is  automated 

in  the  computer  program.  For  delta  wings,  X  =  0,  (x/cr)nt  =  0.667. 

For  a  rectangular  wing,  X  =  1,  equation  (1)  is  used  for  3®  >.  1  and 

the  above  table  is  used  for  3®  <  1.  A  table  of  values  of  (x/c*-) 

r  nt 

for  intermediate  values  of  the  taper  ratio,  X,  obtained  from 
DATCOM  for  the  case  of  no  trailing-edge  sweep  is  included  in  the 
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program.  Linear  interpolation  in  A  is  used  to  obtain  the  value 
of  (x/cr)nt  for  a  particular  fin.  The  values  obtained  from  DATCOM 
are  listed  in  table  1.  The  quantity  tanA  is  the  tangent  of  the 

IjIIj 

leading-edge  sweep  angle  of  the  fin.  If  3/tanATr,  is  less  than 

i-i-L 

1.0,  the  fin  leading  edge  is  subsonic. 

3.3  Axial  Center-of-Pressure  Thickness 
Correction 

The  thickness  distribution  of  a  fin  causes  a  forward  shift 
from  the  linear  theory  value  in  the  axial  center-of-pressure 
position  which,  as  was  shown  in  Reference  1,  can  be  significant. 

For  the  cases  examined  there,  shifts  which  were  as  much  as 
15  percent  of  the  root  chord  were  found.  See  figures  24  through 
33  of  Reference  1.  Thus,  if  an  accurate  prediction  of  fin  hinge 
moment  is  to  be  made,  this  thickness  caused  shift  must  be  accounted 
for.  If  the  leading-edge  shock  wave,  in  a  plane  normal  to  the  fin 
planform  and  parallel  to  the  root  chord,  is  attached,  a  strip- 
theory  method  using  shock-expansion  theory  is  used  to  calculate 
the  center-of-pressure  shift.  No  method  exists  for  accurately 
calculating  the  shift  when  the  leading-edge  shock  wave  is  detached. 
A  method  which  will  yield  a  result,  even  though  it  is  not  valid 
for  this  case,  is  a  strip-theory  method  using  Busemann  second- 
order  theory.  This  method  has  been  used  in  the  present  computer 
program.  The  methods  for  the  attached  and  detached  shock-wave 
cases  will  now  be  described. 

3.3.1  Shock-expansion  theory  for  attached  shock  waves.-  For 
the  attached  shock-wave  case,  strip  theory  is  used  to  calculate 
the  fin  axial  center-of-pressure  shift  due  to  fin  thickness.  In 
this  method  the  fin  is  divided  into  a  series  of  strips  across  the 
span  of  the  fin  and  the  center-of-pressure  shift  using  shock- 
expansion  theory  is  calculated  for  each  strip.  The  shifts  for 
the  various  strips  are  integrated  spanwise  across  the  fin  to 
determine  the  shift  for  the  complete  fin.  The  details  of  the 
method  will  not  be  presented  in  this  report  since  they  are  con¬ 
tained  in  section  4.2.2  and  Appendix  B  of  Reference  1. 
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All  of  the  above  has  been  added  to  the  computer  program. 

The  method  as  programmed  is  restricted  to  the  family  of  airfoils 
described  in  section  3.1  of  this  report.  The  method,  itself,  is 
valid  for  a  fin  planform  composed  of  any  number  of  straight-line 
segments  and  a  varying  airfoil  section  across  the  span  of  the  fin. 
The  computer  program  starts  at  an  angle  of  attack  of  2  degrees 
and  calculates  the  center-of-pressure  shift  due  to  thickness, 
(Ax/cr)t,  up  to  the  angle  of  attack  approximately  equal  to  that 
for  shock  detachment. 


3.3.2  Busemann  theory  for  detached  shock  wave.-  For  the 
detached  shock-wave  case,  strip  theory  is  also  used  to  calculate 
the  fin  axial  center-of-pressure  shift  due  to  fin  thickness.  The 
difference  between  this  case  and  the  attached  shock  case  is  that 
Busemann  second-order  theory  is  used  to  calculate  the  airfoil- 
section  pressure  distribution  rather  than  shock-expansion  theory. 
Busemann  theory  does  not  apply  to  the  detached  shock-wave  case  but, 
as  will  be  seen,  it  does  exhibit  the  correct  behavior  relative  to 
center-of-pressure  values  where  the  shock  wave  is  attached.  In 
lieu  of  a  better  method,  it  is  used  in  the  present  computer  pro¬ 
gram  but  the  results  should  probably  be  used  cautiously. 


The  Busemann  theory  pressure  distribution  is  calculated  using 
the  following  equations  which  are  taken  from  page  243  of  Refer¬ 
ence  6. 


where 


Cp  =  C1 (a  +  ©i)  +  C2(a  +  6i) 

l 


2  ■> 


CP  =  “Cl{a  '  V  +  C2(a  “  V  ^ 
u . 

l 


C,  = 


4 


Moo  -  1 


c„  = 


(Y  +  DM.  -  4(Meo  -  1) 

2(Mj  -  l)2 


(2) 


(3) 
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In  the  above  equations 

a  =  fin  angle  of  attack  (see  fig.  1) 

0jL  =  slope  of  the  i th  region  of  the  airfoil  section  (see 
fig.  1) 

M^  =  free-stream  Mach  number 

Y  =  ratio  of  specific  heats,  y  =  1.4  for  air 

The  subscript  i  refers  to  the  i th  region  of  the  airfoil  section 
as  described  in  Appendix  B  of  Reference  1.  Equations  (2)  and  (3) , 
above,  replace  the  equations  presented  in  section  B.l  of  that 
Appendix  for  calculating  the  airfoil-section  pressure  distribu¬ 
tion.  This  method  is  incorporated  in  the  computer  program  for 
calculating  the  center-of-pressure  shift  due  to  fin  thickness. 

The  calculation  is  only  done  for  a  =  2°  since  the  Busemann  theory 
does  not  apply  to  the  detached  shock  case.  This  value  is  added 

to  the  nonthickness  center  of  pressure,  (x/c  )  ,  and  linear  in- 

r  nt 

terpolation  between  the  resulting  value  and  the  a  =  45°  value, 
obtained  as  described  in  the  following  section,  is  used  to  obtain 
the  center  of  pressure  at  other  angles  of  attack. 

To  show  how  the  center-of-pressure  location  for  a  Mach  number 
where  the  shock  wave  is  detached  fairs  into  the  values  where  it 
is  attached,  calculations  were  made  for  the  fin  shown  in  figure  2. 
The  results  are  shown  in  figure  3.  The  quantity  x/c^  is 


r  —  ^ 

X 

c 

c 

r; 

nt 

r 

v*  y 

For  M  =1.2  and  1.3  the  shock  wave  is  detached  and  for  the  other 

oo 

Mach  numbers  it  is  attached.  Going  from  M^  =  1.6  to  M^  =  1.2 
there  is  a  forward  movement  of  the  center  of  pressure.  If  the 
experimental  data  shown  in  figures  15  through  18  of  Reference  1 
are  examined,  the  same  trend  is  observed.  Thus,  the  Busemann 
theory,  even  though  not  strictly  applicable  to  a  detached  shock 
case,  produces  a  center-of-pressure  shift  in  the  correct  direction. 
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3.4  Axial  Center-of-Pressure  Location 
for  an  Angle  of  Attack  of  45° 

Sections  3.2  and  3.3  present  methods  for  predicting  the  fin- 
alone  center-of-pressure  location  as  a  function  of  angle  of  attack 
up  to  the  'angle  at  which  shock  detachment  occurs.  If  the  shock 
is  detached  at  a  =  2°  a  method  is  presented  for  estimating  the 
center-of-pressure  location  at  this  angle.  Knowledge  of  the  value 
of  x/c^  for  the  fin  at  some  large  angle,  like  45°  where  x/c^  tends 
toward  the  centroid  of  area,  would  allow  the  low  angle  values  to 
be  faired  to  this  value.  In  this  way  an  estimate  of  x/cr  over  a 
large  angle  range  could  be  made. 

To  obtain  values  of  x/c^  at  a  =  45°  for  a  range  of  Mach 
numbers,  M^,  aspect  ratios,  KL,  and  taper  ratios,  A,  the  wing- 
alone  data  base  of  Reference  7  has  been  utilized.  These  data  are 
plotted  in  figure  4  along  with  data  from  the  data  base  of  Refer¬ 
ence  8.  Curves  have  been  faired  through  the  Reference  7  data. 

This  set  was  chosen  since  the  data  of  Reference  8  were  obtained 
on  a  semi-span  model  mounted  on  a  reflection  plane  and  may  be 
affected  by  boundary-layer  separation  at  large  angles  of  attack. 
Values  of  x/c^  have  been  read  from  the  faired  curves  and  are 
tabulated  in  table  2.  These  data  have  been  incorporated  into  the 
computer  program  as  a  data  base.  Triple  linear  interpolation  is 
performed  in  the  data  to  obtain  x/c  at  a  =  45°  for  a  given  Mx, 

and  A.  If  the  value  of  iR  is  greater  than  2.0,  linear  extrapol¬ 
ation  in  A*  is  used. 

3.5  Center-of-Pressure  Calculation 
Procedure 

For  a  given  Mach  number,  the  computer  program  calculates  a 
table  of  x/c^.  versus  a  using  the  methods  described  in  the  previous 
three  sections,  sections  3.2  through  3.4.  The  first  step  in  this 
calculation  is  to  use  the  method  of  section  3.2  to  calculate  the 
fin  alone  center-of-pressure  position  for  the  given  values  of  M  , 
A*,  and  A.  This  provides  the  value  of  (x/cr)  ,  the  center  of 
pressure  without  thickness  effects.  The  thickness  correction, 
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(Ax/c^J^,  is  next  calculated  as  a  function  of  a  for  the  same  value 
of  Mro ,  the  input  airfoil  section,  and  the  fin  planform.  This  cal¬ 
culation  is  made  using  the  procedures  described  in  section  3.3  up 
to  the  shock-wave  detachment  angle.  The  center  of  pressure  of 
the  fin  alone  as  a  function  of  a  is  then  determined  using  Equa¬ 
tion  (4)  .  A  table  of  x/c^.  for  angles  of  attack  up  to  the  leading- 

edge  shock-wave  detachment  angle  has  now  been  determined.  The 
last  entry  in  the  table  is  next  calculated  for  a  =  45°  using  the 
method  in  section  3.4. 

This  table  of  x/cr  versus  a  is  now  used  to  determine  the 

center  of  pressure  of  the  fin  in  the  presence  of  the  body  for  the 

given  values  of  M^,  tf> ,  6,  and  .  The  first  step  in  the  program 
is  to  calculate  the  normal-force  coefficient  of  the  fin  in  the 
presence  of  the  body.  This  is  done  using  the  methods  described 
in  References  2  and  3.  After  this,  the  fin-alone  normal-force 
data  base  in  the  program  is  entered  to  find  the  angle  of  attack 
which  produces  this  normal  force.  Using  this  angle  of  attack, 
the  table  constructed  as  described  in  the  preceeding  paragraph 
is  entered  to  determine  the  value  of  x/c corresponding  to  this 
angle  of  attack.  This  value  of  x/cr  is  used  in  the  hinge-moment 
coefficient  calculation. 

3.6  Maximum  and  Minimum  Hinge 
Moment  Tabulation 

In  order  to  determine  an  optimum  fin  planform  and  airfoil 
section  which  minimizes  fin  hinge  moment,  calculations  must  be 
made  for  a  large  number  of  different  fins  over  the  complete  range 
of  flight  conditions  to  be  encountered,  that  is,  Mach  number, 
angle  of  attack,  and  fin  deflection  angle.  Program  MISSILE2  as 
used  in  References  1  and  3  only  allows  one  Mach  number  and  one 
fin  deflection  angle  per  case  so  that  many  cases  would  have  to  be 
run  to  examine  the  hinge  moment  for  all  flight  conditions  (com¬ 
binations  of  M  ,  4),  a  ,  and  6).  To  alleviate  this,  the  program 
oo  c 

has  been  modified  to  allow  a  series  of  values  of  M  and  6  to  be 

OO 

run  in  the  same  case. 


20 


The  program  has  also  been  changed  so  that  the  maximum  and 
minimum  hinge  moments  are  summarized  in  a  table.  A  sample  of 
this  table  is  shown  in  figure  5.  For  a  given  roll  angle,  the 
maximum  and  minimum  hinge  moments  are  tabulated  as  a  function  of 
Mach  number,  M  ,  and  hinge-line  location  x  /c  ,  for  fin  4. 

ri  IT 

Fin  4  is  the  right  horizontal  fin  viewed  from  the  rear  with  the 
configuration  at  zero  degrees  roll.  The  direction  of  positive 
hinge  moment,  CRM,  is  shown  in  figure  6.  The  values  tabulated 
in  figure  5  are  determined  by  calculating  CHM  for  all  a  ,  6 
combinations  for  a  given  xITT  and  M  .  This  table  of  values  is 
then  searched  to  find  the  maximum  and  minimum  values. 

4.  DETERMINATION  OF  OPTIMUM  HINGE-LINE 
LOCATION  FOR  A  PARTICULAR  FIN 

The  optimum  hinge-line  location  is  the  point  along  the  fin 
root  chord  where  the  magnitudes  of  the  maximum  and  minimum  hinge- 
moment  coefficients  reach  a  minimum  value  when  they  are  considered 
together.  To  illustrate  how  this  is  determined,  the  table  of 
values  in  figure  5  will  be  used.  This  table  was  discussed  in 
section  3.6.  The  body-canard  combination  for  which  the  calcula¬ 
tion  was  made  is  shown  in  figure  7  and  the  fin  details  in  figure  2. 
The  flight  conditions  used  were 

M  =  1.2,  1.6,  2.0,  2.5,  3.0 

CO  r  r  r  f 

4>  =  0° 

61,3  =  °° 

62  4  =  -20°,  -10°,  0°,  10°,  20° 

ac  =  0°,  5°,  10°,  15°,  20° 

This  is  a  total  of  125  points. 

Scanning  down  the  table  we  see  that  both  the  maximum  and 
minimum  hinge  moments  reach  a  minimum  magnitude  at  x  /c  of  0.63 
or  0.64.  The  hinge  moments  in  this  region  will  be  plotted.  To 
determine  the  points  to  plot,  the  maximum  magnitude  for  the  five 
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Mach  numbers  for  which  calculations  were  made  is  selected.  These 
values  are  underlined  in  the  table  from  x„T/c  =  0.60  to  0.67  and 

nJLi  IT 

are  plotted  in  figure  8. 


The  optimum  hinge-line  location  is  the 

If  the 

hinge  line  is  located  at  this  point,  the  maximum  magnitude  of  the 
hinge  moment  for  the  above  flight  conditions  will  not  exceed 


point  at  which  the  two  curves  intersect,  x„T/c  =  0.6322 

nli  r 


'HM 


=  0.015.  Sometimes  the  two  curves  exhibit  the  behavior 


shown  in  the  following  sketch.  In  this  case  the  optimum  location 


of  the  hinge  line  is  not  where  the  two  curves  intersect  but  at 
the  point  where  the  upper  curve  reaches  a  minimum. 

The  above  procedure  has  been  used  in  obtaining  all  of  the 
results  to  be  presented  in  section  5. 

5.  DETERMINATION  OF  OPTIMUM  FIN  TO  MINIMIZE 
HINGE-MOMENT  COEFFICIENT 

The  method  which  has  been  described  will  now  be  used  to 
determine  the  planform  and  the  airfoil  section  of  a  fin  which 
will  minimize  the  magnitude  of  the  fin  hinge-moment  coefficient. 
The  configuration  used  is  the  body  shown  in  figure  7  with  various 
fins.  The  leading  edge  of  the  root  chord  of  all  fins  is  at  the 
location  shown  in  this  figure.  The  body  extends  to  the  fin  trail¬ 
ing  edge.  Since  the  fin  root  chord  will  vary  in  length,  the  body 
length  will  vary  from  that  shown  in  the  figure. 

For  all  of  the  calculations  to  be  shown,  the  flight  condition 
envelope  is  the  same  as  that  used  in  section  4.  That  is 
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1.2,  1.6,  2.0,  2.5,  3.0 


4>  =  o° 

«1,3  -  °° 

52  4  =  -20°,  -10°,  0°,  10°,  20° 

ac  =  0° ,  5° ,  10° ,  15° ,  20° 

The  roll  angle,  <J> ,  is  not  varied  since  experimental  data  show 
that  the  maximum  hinge  moment  acting  on  a  canard  fin  occurs  very 
close  to  <f>  =  0°,  fin  horizontal,  and  is  almost  equal  to  the  value 
at  <f>  =  0°.  This  is  shown  in  figures  84  through  89  of  Reference  1. 
If  strong  body-vortex  effects  are  present  this  may  not  be  the 
case  and  <f>  must  be  varied.  The  same  is  true  when  forward  fins 
induce  large  effects  on  rear  fins.  The  value  of  <J>  can  be  varied 
in  the  computer  program  if  desired. 

The  following  parameters  are  available  for  specifying  the 
fin  planform. 

Aspect  ration,  JR 
Fin  area,  Sp 

Fin  taper  ratio,  X 

Fin  root-chord  length,  c 

Exposed  semispan  of  fin,  (s  -  a) 


In  the  results  to  be  presented,  the  value  of  Sv  will  be  held 

2 

constant  at  14.0625  in  .  Two  aspect  ratios  will  be  considered, 

JR  =  2.0  and  3.0.  The  optimization  will,  thus,  be  done  for  two 
different  aspect  ratio  fins  with  the  same  area.  The  equation  for 
the  aspect  ratio  is 


JR 


4 (s  -  a) 
cr(l  +  X) 


(5) 


where  the  aspect  ratio  is  that  of  two  fins  joined  together  at 
their  root  chords.  The  fin  area  is 


c 

Sp  =  (1  +  X)  (s  -  a) 


(6) 
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From  these  two  equations 


-\/S?  (7 

and 

^SF 

Cr  “  (s  -  a)  (1  +  X)  {8 

Two  parameters  specify  the  airfoil  section 

Thickness  to  chord  ratio,  t/c 
Flat-section  length  to  chord  ratio,  f/c 

For  all  of  the  results,  the  reference  area  is  taken  to  be 

the  cross-sectional  area  of  the  cylindrical  section  of  the  body, 
2 

S  =  19.635  in  ,  and  the  reference  length  is  the  body  diameter, 

£%. 

d  =  5  in.  These  are  used  in  forming  all  force  and  moment 
coefficients . 

The  order  of  optimization  used  is  taper  ratio,  airfoil  sec¬ 
tion,  body-radius  to  fin-span  ratio,  and,  finally,  aspect  ratio. 
These  will  now  be  discussed  in  order. 

5.1  Taper-Ratio  Optimization 

For  the  taper  ratio  optimization  study,  the  fin  airfoil 
section  for  both  JR  =  2.0  and  3.0  was  held  constant  with 

t/c  =0.06 
f/c  =0.25 

The  quantity  (s  -  a)  is  a  function  of  aspect  ratio  and  fin 
area  [Eq.  (7)]  so  that  for  the  two  aspect  ratios 

a/s  =0.4  ;  JR  =  2 . 0 
a/s  =  0.3525  ;  JR  =  3 . 0 

The  value  of  the  magnitude  of  the  maximum  hinge  moment  with 
the  hinge  line  in  the  optimum  position  is  shown  in  figure  9  as  a 
function  of  fin  taper  ratio.  For  both  aspect  ratios,  the  lowest 
hinge  moment  occurs  for  a  fin  of  delta  planform,  X  =  0.  As  a 
result,  a  delta  planform  is  determined  to  provide  the  minimum 
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magnitude  hinge  moment  over  the  envelope  of  flight  conditions. 
A  value  of  A  =  0  will  be  used  in  the  following  steps  of  the 
optimization . 

5.2  Airfoil-Section  Optimization 


The  results  of  the  airfoil  section  optimization  study  are 
shown  in  figure  10  for  the  optimum  taper  ratio  of  A  =  0.  In  this 


figure  the  value  of  J  cHM0pt  Pl°ttec^  against  thickness  ratio, 

t/c,  for  various  values  of  the  flat  plate  parameter,  f/c.  For 
both  aspect  ratios,  the  same  behavior  is  observed.  The  minimum 


value  of 


CHM, 


opt 


is  fairly  insensitive  to  f/c  and,  as  this 
parameter  increases  in  value,  the  minimum  value  occurs  at  smaller 
values  of  t/c.  For  both  aspect  ratios,  the  double-wedge  airfoil, 
f/c  =  0,  is  nearly  the  optimum  one.  As  a  result  of  this  study, 
the  following  airfoil-section  parameters  were  chosen  for  the  two 
aspect  ratios.  They  are  very  near  the  optimum  values. 


m  =  2.0, 

m  =  3.o, 


t/c  =  0.06, 
t/c  =  0.05, 


f/c  =  0.5 
f/c  =0.5 


5.3  Body-Radius  to  Fin-Span  Ratio 
Optimization 

For  a  taper-ratio  zero  fin  and  using  the  above  airfoil- 
section  quantities,  calculations  were  made  varying  the  body-radius 
to  fin-span  ratio,  a/s.  The  results  of  these  calculations  showed 


that  the  value  of 


CHM, 


opt 


did  not  vary  with  a/s . 


5.4  Aspect-Ratio  Optimization 


Results  have  been  presented  in  the  preceeding  sections  for 
values  of  the  aspect  ratio  of  2.0  and  3.0.  Figure  10  shows  that 


the  minimum  value  of 


CHM 


opt 


increases  with  increasing  aspect 
ratio.  This  is  probably  to  be  expected  since  the  normal  force  at 
a  given  angle  of  attack  increases  with  aspect  ratio,  everything 
else  being  held  constant.  This  could  be  verified  by  repeating 
the  above  calculations  for  other  aspect  ratios.  The  normal-force 
coefficients  for  the  two  fins  are  presented  in  the  next  section. 
For  these  minimum  hinge-moment  fins,  control  is  still  maintained. 
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5.5  Results  for  Optimum  Fins 

The  previous  sections  have  taken  two  fins  with  the  same  area, 

2 

S„  =  14.0625  in  ,  and  different  aspect  ratios  and  determined  the 
values  of  the  other  fin  parameters  which  specify  the  planform, 
airfoil  section,  and  hinge-line  location  to  minimize  the  magnitude 
of  the  fin  hinge-moment  coefficient  over  the  range  of  assumed  flight 
conditions.  The  two  fins  selected  are 


Quantity 

Fin  1 

Fin  2 

V  in2 

14.0625 

14.0625 

JR 

2.0 

3.0 

X 

0.0 

o 

• 

o 

t/c 

0.06 

0.05 

f/c 

0.5 

0.5 

a/s 

0.4 

0.3525 

XHl/Cr 

0.6215 

0.6322 

The  computer  program  was  run  in  order  to  examine  the  hinge- 

moment  and  normal-force  coefficient  variation  over  the  range  of 

M  ,  6,  and  a  with  6=0°.  The  results  of  these  calculations  are 
°°  c 

shown  in  figure  11.  Hinge-moment  coefficient,  ,  is  plotted 

versus  normal-force  coefficient,  C.—  .  The  curves  are  for  the 

Nr 

various  free-stream  Mach  numbers. 

The  results  for  the  JR  =  2.0  fin  are  shown  in  figure  11(a). 

The  Ma  =  1.2  curve  has  symbols  showing  the  results  for  the  five 

values  of  a  for  the  five  fin  deflection  angles.  As  can  be  seen 
c 

all  25  points  lie  on  the  same  curve.  Use  of  the  equivalent  angle 
of  attack  concept  causes  this.  This  concept  is  described  in 
section  4.4.4  and  Appendix  C  of  Reference  1,  Appendix  C  of  Refer¬ 
ence  3,  and  Reference  9.  To  obtain  a  desired  CMr,  a  certain  value 

NF 

of  CTT„  is  produced  which  can  be  done  for  various  combinations  of 
HM 

a  and  6.  For  all  combinations  of  a  and  6,  the  normal-force 
c  c 

carryover  onto  the  body  is  the  same.  Since  the  nose  normal 
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force  increases  with  increasing  a  ,  the  a  ,  6  combination  with 

c  c 

the  largest  ac  should  be  used  for  the  best  maneuverability. 

As  can  be  seen,  the  center  of  pressure  for  Mto  =  1.2  is  ahead 

of  the  hinge  line  while  that  for  the  other  values  of  M^  is  behind 
the  hinge  line.  This  behavior  is  consistent  with  experimental 
data  which  show  a  forward  shift  in  the  center  of  pressure  in  the 
transonic  speed  range.  For  this  fin  the  M  =  1.6  to  3.0  results 
are  grouped  quite  close  together. 

The  PR  =  3.0  results  are  shown  in  figure  11(b).  The  results 
are  similar  to  those  for  7R  =  2.0  although  there  is  more  of  a 

Mach  number  effect  for  M  =  1 . 6  to  3 . 0 .  It  is  of  interest  to 

00 

compare  the  results  of  the  two  values  of  /R.  At  M^  =  1.2  the 

maximum  value  of  C>TT1  is  about  the  same  for  the  two  aspect  ratios 

NF 

while  the  maximum  magnitude  of  the  hinge-moment  coefficient  has 
increased  by  about  60  percent  for  the  7R  =  3.0  fin.  At  Mo  =  3.0 
the  maximum  normal-force  coefficient  is  about  140  percent  higher 
for  7R  =  3.0  than  at  2.0  while  the  maximum  magnitude  of  the 
hinge-moment  coefficient  has  increased  by  about  50  percent.  Thus, 
at  M^  =  1.2  for  a  constant  area  fin,  increasing  the  aspect  ratio 
from  2.0  to  3.0  results  in  no  increase  in  normal-force  coefficient 
but  there  is  a  60  percent  increase  in  hinge-moment  coefficient. 

At  =  3.0  the  increase  in  aspect  ratio  produces  a  140  percent 
increase  in  normal-force  coefficient  with  only  a  50  percent 
increase  in  hinge-moment  coefficient. 

6 .  CONCLUDING  REMARKS 

This  report  presents  the  results  of  the  second  phase  of  a 
study  of  the  hinge  moments  of  all-movable  controls  as  used  on 
cruciform  missiles.  In  phase  one  of  the  study,  a  predictive 
method  was  developed  for  predicting  the  hinge  moments  acting  on 
the  controls  and  the  range  of  applicability  of  the  method  was 
determined.  It  was  found  that  the  method  was  applicable  to  the 
Mach  number  range  of  1.3  to  3.7  except  for  leeward  fins  when 
strong  body  vortices  were  present.  As  a  result,  phase  two,  this 
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phase  of  the  study,  was  undertaken  with  the  purpose  of  developing 
an  optimization  technique  using  the  computer  program  which  would 
determine  the  fin  planform,  airfoil  section,  and  hinge-line  loca¬ 
tion  which  would  minimize  fin  hinge-moment  coefficient. 

As  a  result  of  this  goal,  certain  of  the  calculations  done 
in  phase  one  external  to  the  program  have  been  included  in  the 
present  computer  program.  With  minimum  changes  to  the  input  data, 
the  hinge  moments  produced  by  an  all-movable  control  over  a  wide 
range  of  flight  conditions,  fin  planform,  and  fin  airfoil  section 
can  be  studied.  The  configuration  is  restricted  to  a  body  with 
one  set  of  cruciform  fins. 

An  optimization  procedure  is  presented  which  allows  the  user 
of  the  program  to  determine  the  optimum  fin  planform,  airfoil 
section,  and  hinge-line  location  to  minimize  fin  hinge-moment 
coefficient  subject  to  constraints  he  may  apply  through  the  input 
data  to  the  program.  The  range  of  flight  conditions  to  be  studied, 
that  is,  Mach  number,  roll  angle,  angle  of  attack,  and  fin  deflec¬ 
tion  angle,  are  specified  as  input  data.  Through  input  data,  the 
fin  planform  and  airfoil  section  may  also  be  restricted.  In  this 
way  trade-off  studies  can  be  carried  out  over  the  ranges  of 
conditions . 

The  report  presents  an  example  of  an  optimization  study  where 
the  fin  area  is  fixed  and  all  other  parameters  are  free  to  vary 
for  a  given  range  of  flight  conditions.  This  example  illustrates 
the  use  of  the  program  for  a  study  such  as  this.  The  computer 
program  and  its  use  is  described  in  an  appendix  of  this  report. 

The  range  of  applicability  of  the  program  is 

Mach  number;  1.2  to  4.0 
Roll  angle;  -90°  to  90° 

Angle  of  attack;  0°  to  20° 

Fin  deflection  angle;  -20°  to  20° 

Aspect  ratio;  0.5  to  4.0 
Taper  ratio;  0.0  to  1.0 
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The  fin  airfoil  section  is  restricted  to  a  flat  plate  with  equal 
length  wedges  for  the  leading  and  trailing  edges.  This  is  not  a 
restriction  on  the  method  but  one  imposed  by  the  program. 

As  a  result  of  the  optimization  study  presented  in  this 
report,  the  following  conclusions  can  be  drawn. 


1.  Hinge-line  location  was  the  most  important  parameter  in 


minimizing 


CHM, 


opt 


Fin  taper  ratio  had  a  large  effect  with  the  smallest 


CHM, 


occurring  at  low  values  of  X  (0  to  0.25) 


opt 

The  value  of  Cjjm 
fin-span  ratio 

The  minimum  value  of 


opt 


was  insensitive  to  body-radius  to 


CHM, 


opt 


was  quite  insensitive  to 


various  combinations  of  t/c  and  f /c .  However,  if  a 
given  value  of  t/c  must  be  used  for  structural  reasons, 
the  right  value  of  f/c  must  be  selected. 

Increasing  the  aspect  ratio  from  2.0  to  3.0  caused  an 


increase  in  the  minimum 
At  M„  =  1.2  the  maximum  C. 


opt 


of  50  to  60  percent. 


NF 


was  not  changed  whereas  at 


M 


=  3.0  it  was  increased  about  140  percent. 


The  present  work  used  one  optimization  scenario,  that  of 
minimizing  the  maximum  value  of  C„..  for  a  constant  area  fin. 

nM 

The  computer  program  can  be  easily  changed  to  handle  other 
scenarios . 
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Table  1.-  Values  of  (x/cr)nt  obtained  from  DATCOM  for 
>  1.0  for  fins  with  unswept  trailing  edges. 


tanA 

LE 

6 

x/c 

nt 

X  =  .2 

A  =  .25 

X  =  .33 

to 

ll 

o 

o 

0.650 

0.652 

0.640 

0.612 

CM 

o 

0.648 

0.652 

0.635 

0.601 

0.4 

0.645 

0.650 

0.627 

0.585 

0.6 

0.640 

0.644 

0.618 

0.563 

CO 

o 

0.635 

0.630 

0.600 

0.538 

1.0 

0.630 

0.605 

0.578 

0.510 

e 

f  >| 

x/c 

rJ 

nt 

tanA 

LE 

X  =  .2 

X  =  .25 

X  =  .33 

X  =  .5 

o 

o 

0.535 

0.500 

0.450 

0.345 

0.2 

0.575 

0.535 

0.500 

0.386 

0.4 

0.605 

0.565 

0.535 

0.422 

0.6 

0.618 

0.590 

0.557 

0.454 

00 

• 

o 

0.625 

0.600 

0.572 

0.483 

1.0 

0.630 

0.605 

0.578 

0.510 

Table  2.-  Values  of  x/cr  at  a  =  45°  as  determined 
from  the  Stallings-Lamb  data  base. 


M 

CO 

IR  =  1/2 

JR  =  1 

m  =  2 

X  =  0 

X  =  1/2 

X  =  1 

X  =  0 

X  =  1/2 

X  =  1 

X  =  o 

X  =  1/2 

X  =  1 

1.0 

0.601 

0.535 

0.435 

0.570 

0.545 

0.430 

0.618 

0.545 

0.417 

1.4 

0 .626 

0.551 

0.435 

0.610 

0.555 

0.430 

0.620 

0.546 

0.417 

1.8 

0.638 

0.562 

0.436 

0.629 

0.560 

0.430 

0.621 

0.547 

0.416 

2.2 

0.642 

0.570 

0.437 

0.637 

0.565 

0.429 

0.623 

0.548 

0.415 

2.6 

0.643 

0.574 

0.438 

0.640 

0.565 

0.428 

0.625 

0.549 

0.415 

3.0 

0.644 

0.578 

0.439 

0.640 

0.565 

0.428 

0.627 

0.549 

0.414 

3.8 

0.645 

0.580 

0.440 

0.640 

0.565 

0.427 

0.630 

0.550 

0.413 

4.6 

0.646 

0.580 

0.441 

0.640 

0.565 

0.426 

0.633 

0.551 

0.411 

OJ 

U) 


Figure  1.-  Airfoil  section  in  computer  program. 
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Figure  2.-  Dimensions  of  an  £R=3.0,  A=0.0  fin. 
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Figure  3.-  Center-of-pressure  location  as  a  function  of  Mach  number 

for  an  PR  =  3.0,  X  =  0.0  fin. 


(a)  PR  =  1/2 


Figure  4.-  Comparison  of  axial  center-of-pressure  position  at  a  =  45°  of 

the  data  bases  of  references  7  and  8 . 


(b)  m  =  i 


Figure  4.-  Continued. 
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Figure  4.-  Concluded. 


maximum  and  MINIMUM  HINGE  MOMENTS  FOR  F 1 N  4  (BASED  ON  SROU!  AND  LROUT).  UVE«  THE  RANGES  OF  ANGLE  OF  ATTACK 
AND  FIN  DEFLFCTION  ANGLE  FO^  WHICH  CALCULATIONS  WFRE  MADE*  TUB  A  GIVEN  BOLL  ANGLE.  THEY  ARE  TABULATED 

as  a  function  of  hinge-line  location  and  mach  number. 
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Figure  6.-  Positive  directions  for  canard  fin  force 
and  moment  coefficients. 
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Figure  7.-  Body-canard  combination. 


Figure  8.-  Determination  of  optimum  hinge-line  location. 
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(a)  PR  =  2.0 

Figure' 9.-  Variation  with  taper  ratio  of  the  magnitude  of 
the  maximum  hinge-moment  coefficient  with 
the  hinge  line  in  the  optimum  location. 
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(b)  m  =  3.0 


Figure  9.-  Concluded. 
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t/c 

(a)  m  =  2.0 


Figure  10.-  Variation  with  airfoil  section  parameters  of  the 
magnitude  of  the  maximum  hinge-moment  coefficient  with  the 
hinge  line  in  the  optimum  location. 


Figure  10,-  Concluded 


(a)  PR  =  2.0 


Figure  11.-  Variation  of  hinge -moment  coefficient  with 
fin  normal-force  coefficient  for  the  optimum 

fins,  <j>  =  0°. 
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Figure  11.-  Concluded. 
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APPENDIX  A 


USER'S  MANUAL  FOR  COMPUTER  PROGRAM 

A . 1  INTRODUCTION 

The  purpose  of  this  appendix  is  to  describe  the  operation  of 
the  computer  code  in  sufficient  detail  to  permit  understanding 
and  use  of  the  program.  The  program  computes  forces  and  moments 
for  each  section  of  the  configuration,  and  for  the  complete  con¬ 
figuration.  Minimum  drag  is  not  computed.  The  code  is  capable 
of  handling  a  body-canard  or  body-tail  with  no  afterbody.  A 
description  of  the  engineering  method  is  given  in  References  Al , 

A2 ,  and  A3  and  in  the  main  text  of  this  report. 

The  program  is  written  in  FORTRAN  IV  and  has  been  run  on  the 
CDC  7600  and  Cyber  760  machines.  A  typical  running  time  for  five 
Mach  numbers,  five  canard  deflection  angles,  five  angles  of  attack, 
and  one  roll  angle  (125  points)  is  about  13  seconds  on  the  CDC 
7600  and  about  19  seconds  on  the  Cyber  760. 

A. 2  PROGRAM  DESCRIPTION 

For  computational  purposes,  the  missile  to  be  studied  is 
divided  into  two  sections.  The  nose  section  is  defined  to  be 
from  the  nose  tip  to  the  leading  edge  of  the  root  chord  of  the 
finned  section.  The  canard  section  (finned  section)  is  defined 
to  be  from  the  end  of  the  nose  section  to  the  trailing  edge  of 
the  set  of  fins.  The  main  program  controls  the  flow  of  opera¬ 
tions.  It  calls  the  routines  which  control  the  operations  for 
each  section  of  the  missile.  Figure  A.l  shows  the  routines 
called  by  the  main  program  and  the  routines  which  they  in  turn 
call . 
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A. 2.1  Calculation  Procedure 


Figure  A. 2  shows  the  flow  of  the  program  computation.  This 
figure  and  figure  A.l,  which  shows  the  subprogram  call  sequence, 
are  the  basis  for  the  following  description  of  the  calculation 
procedure. 

Subroutine  INPT  reads  in  the  run  identification  information, 
and  then  reads  in  body  and  fin  geometry,  run  control  parameters, 
and  flow  conditions.  All  length  quantities  are  then  normalized 
by  the  body  radius  and  all  areas  are  normalized  by  the  square  of 
the  body  radius . 

After  all  the  input  has  been  read  in,  the  required  inter¬ 
ference  factors  are  computed  by  calling  routines  INTFAC,  INFLU, 
and  CCL .  Next,  routine  WNGCNW  is  called  to  obtain  the  wing-alone 
normal-force  coefficient,  CN ,  as  a  function  of  angle  of  attack, 
a,  from  the  data  base,  for  the  input  canard  planform.  Subroutines 
SHKEXP  and  XCPVSA  are  called  to  calculate  the  wing-alone  longi¬ 
tudinal  center  of  pressure,  x/c ,  as  a  function  of  a.  This 
calculation  includes  the  correction  for  thickness  effects. 

The  next  step  in  the  calculation  is  to  calculate  the  loads 
on  the  nose  section.  This  is  accomplished  through  calls  to  LNTRP 
and  NOSE.  LNTRP  is  used  to  determine  some  crossflow-drag  quanti¬ 
ties.  NOSE  first  computes  the  axial  starting  location  of  the 
nose  vortices  (if  present) .  The  locations  and  strengths  of  the 
nose  vortices  in  the  crossflow  plane  at  the  leading  edge  of  the 
root  chord  of  the  canard  section  are  computed  next.  The  subrou¬ 
tine  then  computes  the  forces  and  moments  on  the  nose  and  returns 
control  to  the  main  program.  The  vortex  positions  and  the  forces 
and  moments  are  printed  in  NOSE  if  the  parameter  NOUTPT  =  1 . 

Subroutine  CANARD  is  now  called  to  calculate  the  canard 
loads  and  the  body  loads  in  the  canard  section.  The  individual 
fin  forces  and  moments  in  the  absence  of  the  nose  vortices  are 
computed  first.  The  user  may  choose  to  ignore  the  nose  vortices 
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over  the  canards  if  he  believes  they  have  dispersed.  This  con¬ 
trol  is  achieved  throught  the  parameter  NVORT. 

If  nose  vortices  are  present  over  the  finned  section,  sub¬ 
routine  REVFLO  is  then  called  by  CANARD.  Subroutine  REVFLO 
computes  the  equivalent  angle  of  attack  and  spanwise  location  of 
the  center  of  pressure  for  the  loading  due  to  the  nose  vortices. 
Subroutine  REVFLO  assumes  the  vortices  to  be  infinite  line 
vortices  parallel  to  the  body  axis. 

When  REVFLO  returns  control  to  CANARD,  the  total  fin  forces 
and  moments  and  body  forces  and  moments  (for  the  canard  section) 
are  computed.  The  canard  individual  fin  loads  are  always  printed 
in  CANARD.  If  the  parameter  NOUTPT  =  1,  additional  output  is 
printed.  This  includes  the  effect  of  hinge-line  location  on  the 
fin  hinge  moments  and  the  contribution  of  the  canard  section  to 
the  total  loads. 

Upon  returning  to  the  main  program,  a  summary  of  the  overall 
forces  and  moments  is  printed  if  NOUTPT  =  1.  Following  this,  a 
search  through  all  canard  deflection  angle,  6,  and  angle  of 
attack,  ac,  combinations  for  this  Mach  number,  Mo ,  and  roll  angle, 
4> /  is  made  to  find,  as  a  function  of  hinge-line  location,  the 
maximum  and  minimum  hinge  moments  acting  on  fin  4.  This  fin  is 
the  right-horizontal  fin,  looking  upstream  at  <J>  =  0°.  When  $  =  0° 
the  fins  are  vertical  and  horizontal. 

After  exiting  from  the  Mach  number  loop,  the  maximum  and 
minimum  hinge  moments  acting  on  fin  4  are  tabulated  as  a  function 
of  hinge-line  location  and  Mach  number.  A  separate  table  is 
printed  for  each  roll  angle. 

A. 2. 2  Program  Limitations  and  Precautions 

The  program  makes  a  number  of  assumptions  about  the  missile 
configuration  and  the  flow  field.  These  are  described  briefly 
below. 
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1.  The  fins  in  the  canard  section  must  be  identical,  uncam¬ 
bered  and  untwisted.  Thickness  effects  are  included  and  the  fin 
airfoil  sections  must  be  similar  across  the  span  of  the  fin. 

2.  The  leading  edges  of  the  fins  must  not  be  swept  forward 
and  the  trailing  edges  must  be  unswept.  If  a  fin  with  rectangular 
planform  is  to  be  modeled  the  user  must  set  the  leading-edge  sweep 
to  some  small  positive  value.  The  procedure  for  doing  this  is 
described  in  section  A. 3. 

3.  The  included  angle  of  attack,  ac,  should  not  exceed  45°. 

4.  The  fin  equivalent  angle  of  attack  should  not  exceed  60°. 

5.  The  Mach  number  range  of  the  data  base  in  the  program 
for  fin  normal-force  coefficient  is  0.8  to  3.0.  The  program  can¬ 
not  be  run  below  a  Mach  number  of  0.8  and  probably  should  not  be 
run  below  M^  =  1.2  because  of  strong  transonic  effects  not 
accounted  for  in  the  methods.  In  the  work  of  reference  A3,  the 
Mach  number  extrapolation  scheme  was  changed  so  the  extrapolation 


If  was  found  there  that,  using  data  in 


the  Ma  =  2  to  3  range,  the  normal-force  coefficient  at  M^  =  4.6 
could  be  well  predicted.  Thus,  the  present  program  can  probably 
be  used  up  to  M  =4.5  with  some  confidence. 

6.  The  aspect  ratio  range  of  the  data  base  in  0.5  to  3.53 
for  0.8  <.  M^  <.  1.3  and  0.5  to  2.0  for  1.3  <_  M^  <,  3.0.  The  program 
will  extrapolate  beyond  the  data  base  to  include  fins  with  aspect 
ratio  less  than  0.5  but  it  is  recommended  that  such  fins  not  be 
used.  This  is  because  of  the  large  variation  in  loads  with  aspect 
ratio  in  this  range.  In  the  work  of  Reference  A3  the  aspect  ratio 
extrapolation  to  aspect  ratios,  JR,  greater  than  2.0  was  changed 
from  linearly  extrapolating  in  JR  to  extrapolating  as  1/JR.  It 
was  found  there  that  the  normal-force  coefficient  for  an  JR  =  4.0 
fin  could  be  predicted  quite  well  using  JR  =1.0  and  2.0  data. 
Therefore,  the  present  program  can  probably  be  used  for  aspect 
ratios  up  to  4.0. 
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A. 2. 3  Description  of  Subroutines 


The  main  program  primarily  acts  as  an  executive  routine. 

It  calls  the  other  subroutines  as  needed  and  totals  the  loads 
from  each  section. 

Subroutine  BUSEMN  uses  Busemann ' s  second-order  theory  to 
calculate  the  center-of-pressure  shift  due  to  the  fin  airfoil 
thickness  distribution.  It  is  described  on  page  243  of 
Reference  A4 . 

Subroutine  BVTEX  computes  nose  vortex  positions  and  strengths 

Subroutine  CANARD  controls  the  subroutines  which  compute 
vortex  strengths  and  positions,  individual  fin  forces  and  moments, 
body  forces  and  moments,  and  total  forces  and  moments  for  the 
canard  section. 

Subroutine  CCL  computes  the  fin  influence  coefficients  used 
in  the  reverse  flow  procedure. 

Function  CHRT8  calculates  the  fin  lift-curve  slope  at  super¬ 
sonic  speeds  from  the  curves  of  Chart  8  of  Reference  A5 .  This 
subroutine,  as  well  as  some  of  the  other  subroutines  are  docu¬ 
mented  in  Reference  A6 . 

Subroutine  CH1416  calculates  the  center  of  pressure  of  the 
lift  carryover  onto  the  body  due  to  the  canard  from  the  curves  of 
Charts  14,  15,  and  16  of  Reference  A5 .  This  subroutine  is 
documented  in  Reference  A6 . 

Subroutine  CLAM  calculates  the  fin-fin  interference  factors 
used  in  determining  the  change  in  equivalent  angle  of  attack  due 
to  fin  deflection. 

Functions  CNT6 ,  CNTll,  CNT14 ,  CNT15 ,  CNT23 ,  and  CNT31  compute 
the  vortex-free  normal-force  coefficient  for  fins  in  the  data  base 
These  data  are  contained  in  tables  8  and  9  of  Reference  Al. 

Subroutine  CURVES  is  a  block  data  routine  for  initializing  a 
number  of  empirical  tables  used  in  the  program. 
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Function  EQ30  calculates  the  value  of  Kg ^ ( 3CNa) c  using 
equation  (30)  of  Reference  A5 .  This  function  is  used  for  the 
high-aspect-ratio  range  at  supersonic  speeds  when  there  is  no 
afterbody  behind  the  finned  section  and  3m  is  greater  than  one. 
This  subroutine  is  documented  in  Reference  A6 . 

Function  EQ31  is  similar  to  EQ30  and  is  used  when  3m  is 
less  than  one. 

Subroutines  FINTML,  FINTMR ,  FINTNL ,  and  FINTNR  compute 
integrands  used  in  the  reverse-flow  procedure  by  CCL.  These 
integrands  are  discussed  in  Appendices  A  and  B  of  Reference  Al . 

Subroutines  IBCIEU,  ICSEVU,  and  ICSICU  are  routines  used 
for  cubic  spline  interpolation  in  the  data  base. 

Subroutine  INFLU  computes  the  effects  of  Mach  number  on 
panel-panel  interference. 

Subroutine  INPT  reads  and  prints  all  input  data  and  non- 
dimensionalizes  it. 

Subroutine  INT  performs  linear  interpolation  in  a  three 
dimensional  array  of  data. 

Subroutine  INTFAC  calculates  interference  factors  by  the 
method  of  slender  body  theory  or  the  methods  of  Reference  A5 . 

Subroutine  LNTRP  is  a  linear  interpolation  routine. 

Subroutine  NOSE  computes  forces  and  moments  on  the  nose  sec¬ 
tion  and  the  positions  and  strengths  of  the  nose  vortices,  if 
any,  at  the  canard  root-chord  leading  edge. 

Subroutine  REVFLO  computes  the  equivalent  angles  of  attack 
and  spanwise  locations  of  the  centers  of  pressure  for  a  set  of 
fins  due  to  the  presence  of  vortices. 

Subroutine  SHKAGL  calculates  the  shock-wave  angle  associated 
with  a  wedge  of  a  given  angle  using  the  method  of  Appendix  B  of 
Reference  A3. 
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Subroutine  SHKEXP  calculates  the  center-of-pressure  shift 
due  to  fin  thickness  using  shock-expansion  theory  by  the  method 
of  Appendix  B  of  Reference  A3. 

Subroutine  SIMPl  is  a  Simpson's  Rule  integration  package 
used  by  CCL. 

Subroutine  SIMSON  is  a  Simpson's  Rule  integration  package 
used  by  REVFLO . 

Subroutine  VEL  calculates  the  velocity  at  several  points 
spanwise  along  a  fin  induced  by  external  vortices. 

Function  WNGCNT  computes  the  vortex-free,  normal-force 
coefficient  for  a  general  fin  from  the  data  base  contained  in 
functions  CNT6 ,  CNTll,  CNT14 ,  CNT15 ,  CNT23 ,  and  CNT31 . 

Subroutine  WNGCNW  computes  the  wing-alone  normal-force  coef¬ 
ficient  for  a  fin.  The  data  used  are  contained  in  tables  2  and  3 
of  Reference  Al . 

Function  XBAR  computes  the  chordwise  location  of  the  fin 
center  of  pressure  from  the  table  calculated  in  subroutine  XCPVSA. 

Subroutine  XCPVSA  calculates  a  table  of  chordwise  locations 
of  the  fin  center  of  pressure  as  a  function  of  angle  of  attack 
including  thickness  effects. 

Function  YTAB  computes  the  spanwise  location  of  the  fin 
center  of  pressure  from  the  data  base  given  in  tables  4  and  5  of 
Reference  Al . 

A  chart  listing  which  common  blocks  are  in  each  subroutine 
is  given  in  figure  A. 3. 

A. 3  DESCRIPTION  OF  INPUT 

This  section  describes  the  input  required  by  the  computer 
program.  Included  is  a  discussion  of  any  constraints  on  the  value 
or  use  of  the  variables  and,  where  appropriate,  suggested  values. 
Basically,  there  are  three  types  of  variables;  geometric,  flow. 
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and  program  control.  In  addition,  the  order  in  which  the  input 
variables  are  read  approximately  corresponds  to  the  geometric 
order  of  the  various  sections  of  a  missile.  That  is,  the  infor¬ 
mation  required  by  the  nose  section  is  entered  first,  followed  by 
the  information  for  the  canard  section. 

The  program  has  been  designed  to  study  several  different 
configurations  consecutively.  Thus,  the  input  decks  for  each  of 
the  configurations  can  be  stacked  together,  with  the  exception 
of  item  12,  and  the  program  will  analyze  each  configuration  in 
order.  Item  12  is  the  last  card  of  the  input  deck  and  indicates 
the  end  of  information.  All  input  variables  are  listed  and  de¬ 
fined  in  the  next  section  in  the  order  of  appearance  in  the  input 
deck.  The  input  formats  are  shown  in  figure  A. 4.  The  item  num¬ 
bers  below  also  refer  to  that  figure.  A  sample  input  case  is 
discussed  in  section  A. 6. 

Item  1 

These  cards  provide  identification  of  the  run.  The  first 
card  contains  the  index  NCARDS  which  indicates  how  many  cards  of 
information  follow  to  describe  the  run.  The  value  of  NCARDS  must 
be  one  or  greater. 

Item  2 

This  item  provides  run  control  information.  The  first  vari¬ 
able,  NMACH ,  specifies  the  number  of  values  of  the  Mach  number, 

M  ,  for  which  calculations  are  to  be  made.  The  value  of  NMACH 
must  be  between  one  and  five. 

The  second  variable,  NDELT ,  specifies  the  number  of  sets  of 
canard  fin  deflection  angles,  6^,  i  =  1,2, 3, 4,  for  which  calcula¬ 
tions  are  to  be  made.  The  value  of  NDELT  must  be  between  one  and 
five . 

The  third  variable,  NALFA,  specifies  the  number  of  values  of 
the  angle  of  attack,  ac ,  at  which  calculations  are  to  be  made. 

The  value  of  NALFA  must  be  between  one  and  ten. 
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The  fourth  variable,  NPHI,  indicates  the  number  of  roll 
angles,  4> ,  for  which  calculations  are  to  be  made.  The  value  of 
NPHI  must  be  between  one  and  ten. 

The  computer  program  consists  of  a  quadruple  do  loop  as 
shown  in  figure  A. 2.  Thus,  the  total  number  of  conditions  for 
which  calculations  are  made  is  the  product  of  these  four 
variables . 

The  fifth  variable,  NI ,  is  the  number  of  integration  inter¬ 
vals  plus  one,  used  in  the  Simpson's  rule  integration  routine 
called  by  REVFLO.  The  value  of  NI  must  be  odd.  A  suitable  value 
for  NI  depends  on  how  close  a  vortex  is  to  a  fin.  If  a  vortex 
is  close  to  a  fin,  the  value  of  NI  should  be  at  lease  51.  Since 
REVFLO  uses  only  a  small  fraction  of  the  total  run  time,  it  is 
recommended  that  NI  be  set  equal  to  99,  the  maximum  possible 
value . 

The  sixth  variable,  NNOSE ,  is  equal  to  the  number  of  entries 
in  the  table  of  nose  coordinates.  The  first  value  is  the  dis¬ 
tance  of  the  nose  tip  from  missile  station  zero  and  the  last 
value  is  the  location  of  the  shoulder.  The  value  of  NNOSE  must 
be  between  two  and  twenty. 

The  seventh  variable,  NCA,  indicates  whether  the  slope  of 
the  linear  normal-force  curve,  Cfla,  for  the  nose  is  to  be  read 
in.  If  NCA  equals  zero,  then  CNa  is  assigned  the  value  2.0  by 
the  program.  This  is  the  value  predicted  by  slender-body  theory. 
If  NCA  ^  0,  a  value  for  CNa  is  to  be  read  in  later. 

The  eighth  variable,  NVORT,  is  used  to  control  the  influence 
of  the  nose  vortices.  It  has  been  observed  (Refs.  A7  and  A8) 
that  nose  vorticity  for  ac  £  20°  may  disperse  over  the  canard 
section.  Since  the  present  model  for  nose  vorticity  is  incapable 
of  representing  such  a  situation,  the  user  has  the  option  of 
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ignoring  the  influence  of  nose  vorticity  downstream  of  the  canard 
root-chord  leading  edge.  The  options  are: 

NVORT  =  1:  nose  vortices  are  ignored  downstream  of  the 
leading  edge  of  the  canard  root  chord 

NVORT  =  2:  nose  vortices  are  included  to  the  trailing  edge 
of  the  canard 

The  ninth  variable,  NOUTPT,  controls  the  amount  of  output 
produced  by  the  program.  If  NOUTPT  =  0,  the  canard  fin  forces  and 
moments  are  tabulated  as  well  as  the  final  summary  tables  of  the 
maximum  and  minimum  hinge  moments.  If  NOUTPT  =  1,  the  above  is 
printed  as  well  as  the  nose  vortex  strengths  and  positions,  the 
nose  section  forces  and  moments,  additional  canard  fin  informa¬ 
tion,  the  total  canard  section  forces  and  moments,  and  the  total 
configuration  forces  and  moments. 

The  tenth  variable,  logical  variable  TURB,  is  used  to  deter¬ 
mine  which  branch  of  the  crossflow  drag  coefficient  table  is  to 
be  used.  For  crossflow  Mach  numbers  below  0.6,  the  laminar  and 
turbulent  values  of  c^  differ.  If  there  is  doubt  about  which 
type  of  flow  separation  is  present,  use  TURB=.TRUE.,  turbulent 
flow. 

The  eleventh  variable,  logical  variable  REFER,  is  used  to 
control  the  reference  areas  and  lengths  for  the  fin  output.  If 
REFER  =  .TRUE.,  the  reference  area  for  the  fin  loads  is  the  fin 
planform  area,  the  reference  length  for  fin  hinge  moments  is  the 
fin  root  chord,  and  the  reference  length  for  fin  root-bending 
moments  and  rolling  moments  is  the  exposed  span  of  the  fin.  If 
REFER  =  .FALSE.,  the  input  reference  area  and  length  are  used  for 
the  fin  output  as  well  as  for  the  overall  loads. 

Item  3 

This  item  provides  some  geometric  information.  The  first 
two  variables,  LROUT  and  SROUT,  are  the  reference  length  and 
reference  area,  respectively.  The  next  variable  is  XMC ,  the 
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moment  center  measured  from  missile  station  zero  (MSO) .  The 
variable  A  is  the  radius  of  the  cylindrical  section  of  the 
missile . 

The  variables  TIPRAD  and  ETAN  are  entered  next.  The  vari¬ 
able  TIPRAD  is  the  radius  of  the  spherical  nose  tip.  If  the  nose 
is  pointed,  then  ETAN  is  the  nose  half-angle.  If  the  nose  is 
blunted,  then  ETAN  is  the  angle  between  the  body  axis  and  the 
tangent  to  the  nose  at  the  juncture  of  the  spherical  cap  and  the 
rest  of  the  nose. 

The  variable  XCLE  is  the  distance  from  missile  station  zero 
(MSO)  to  the  leading  edge  of  the  root  chord  of  the  canard  fins. 
Any  dimensional  system  is  acceptable.  However,  one  must  be  care¬ 
ful  to  use  the  chosen  system  consistently. 

Item  4 

This  item  contains  the  NMACH  values  of  the  Mach  number  for 
which  computations  are  to  be  performed. 

Item  5 

The  variable  DCNDA  is  the  slope  of  the  nose  normal-force 
coefficient  curve  due  to  attached  flow.  The  reference  area  is 
the  base  of  the  nose.  This  item  is  included  only  if  the  value 
of  NCA  in  Item  2  was  set  equal  to  a  value  different  than  zero. 

Item  6 

This  item  contains  the  nose  coordinates.  First  the  axial 
positions,  XNOSE ,  are  entered,  followed  by  the  corresponding 
values  of  the  local  nose  radius,  RNOSE .  There  should  be  NNOSE 
values  of  each.  The  values  of  XNOSE  are  measured  from  missile 
station  zero. 

Item  7 

This  item  contains  the  NALFA  values  of  the  angle  of  attack 
for  which  computations  are  to  be  performed. 
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Item  8 


This  item  contains  the  NPHI  values  of  the  role  angle  for 
which  computations  are  to  be  made. 

Item  9 

This  item  is  a  set  of  NDELT  cards.  Each  card  of  this  item 
contains  the  fin  deflection  angles  for  the  four  canard  fins.  As 
viewed  from  the  rear,  for  <j>  =  0°,  6^,  and  6^  (vertical  fins)  are 
positive  for  trailing  edges  to  the  right  and  62  and  <5^  (horizontal 
fins)  are  positive  for  trailing  edges  down.  At  <p  =  0°,  fin  1  is 
on  top  and  the  fins  are  numbered  counterclockwise. 

Item  10 

This  item  is  one  card  and  contains  geometric  information 
about  the  canard  section.  The  first  variable,  SPANC,  is  the  fin 
semispan  measured  from  the  body  axis.  The  next  three  variables, 
XCHL,  XCTIPL,  and  XCTE ,  are  the  axial  positions  of  the  hinge  line, 
the  leading  edge  of  the  tip  chord,  and  the  trailing  edge,  respec¬ 
tively.  These  three  axial  positions,  as  well  as  all  other  axial 
positions,  are  measured  from  MSO.  Recall  from  section  A. 2. 2  that 
the  leading-edge  sweep  cannot  be  zero,  even  for  a  rectangular 
planform.  For  the  case  of  a  rectangular  planform  set  XCTIPL 
greater  than  XCLE  by  a  small  number,  for  example  0.01. 

Item  11 

This  item  contains  the  data  required  to  calculate  the  chord- 
wise  center-of-pressure  shift  due  to  fin  thickness.  The  first 
variable,  NFOILS,  is  the  number  of  strips  across  the  span  which 
the  fin  is  to  be  divided  into.  It  must  be  between  one  and  ten. 

For  a  rectangular  fin,  one  can  be  used.  For  other  fins,  five  or 
six  is  sufficient.  The  last  two  variables  of  this  item,  TOC  and 
FOC,  specify  the  airfoil  section  which  is  similar  across  the  span 
of  the  fin.  The  program  treats  an  airfoil  which  is  a  flat  plate 
with  equal  length  leading  and  trailing  edge  wedges,  the  variable 
TOC  is  the  ratio  of  the  flat-plate  thickness  to  the  airfoil  chord 
while  FOC  is  the  ratio  of  the  length  of  the  flat  section  to  the 
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airfoil  chord.  FOC  can  be  zero,  a  double  wedge  airfoil,  and  can¬ 
not  be  one . 


This  is  the  last  card  of  input  for  a  case.  If  an  additional 
case(s)  is  to  be  included  in  this  run,  start  over  with  Item  1. 
Otherwise  add  Item  12. 

Item  12 

This  card  ends  the  process  of  entering  data.  It  has  999 
punched  in  columns  3,  4,  and  5.  It  should  be  the  last  card  and 
follow  all  the  data  cards  for  the  case(s)  to  be  run.  The  com¬ 
puter  program  stops  the  search  for  more  data  and  the  run  is 
completed . 


A. 4  SYMBOL  LISTING 

Program 

Variable 

Item  1 

Alphanumeric  information  to  identify  the  run. 

NCARDS 

Number  of  cards  used  to  identify  the  run; 

NCARDS  >_  1. 

HEAD (I) 

NCARDS  cards  of  alphanumeric  information  for  iden¬ 
tification  of  the  run;  1  <  I  <  NCARDS. 

Item  2 

Integer  and  logical  variables  for  control  of 
program  operation. 

NMACH 

Number  of  Mach  numbers  for  which  calculations  are 
to  be  made ;  1  <  NMACH  <  5 . 

NDELT 

Number  of  sets  of  canard  fin  deflection  angles  for 
which  calculations  are  to  be  made;  1  <  NDELT  £  5. 

NALFA 

Number  of  angles  of  attack  for  which  calculations 
are  to  be  made;  1  <  NALFA  <  10. 

NPHI 

Number  of  roll  angles  for  which  calculations  are 
to  be  made;  1  <  NPHI  <  10. 

NI 

One  plus  the  number  of  intervals  to  be  used  in  the 
Simpson's  rule  integration  package  in  REVFLO;  must 
be  odd;  1  <  NI  <  99. 
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Program 

Variable 


NNOSE 

Number  of  entries  in  the  table  of  nose  coordinates; 

2  <  NNOSE  £  20. 

NCA 

Integer  flag  specifying  whether  Cjga  of  nose  is  to 
be  entered . 

NCA  =0;  Cm  not  entered. 
lNa 

NCA  ^  0:  CNa  entered. 

NVORT 

Integer  flag  indicating  how  far  along  body  influence 
of  nose  vortices  is  to  be  felt. 

NVORT  =  1:  influence  of  nose  vortices  felt  up 
to  leading  edge  of  canard  root 
chord . 

NVORT  =  2:  influence  of  nose  vortices  felt  up 
to  trailing  edge  of  canard  root 
chord . 

NOUTPT 

Integer  flag  controlling  amount  of  output. 

NOUTPT  =  0:  only  fin  loads  and  maximum  and 

minimum  hinge  moments  are  output. 
NOUTPT  =  1:  complete  output. 

TURB 

Logical  variable  stating  whether  crossflow  on  body 
is  laminar  or  turbulent. 

TURB  =  .TRUE.:  crossflow  is  turbulent. 

TURB  =  .FALSE.:  crossflow  is  laminar. 

REFER 

Logical  variable  concerning  output  reference  areas 
and  lengths  for  fins. 

REFER  =  .TRUE.:  use  fin  planform  area,  root 

chord  and  exposed  semispan. 
REFER  =  .FALSE.:  use  input  reference  area  and 

reference  length. 

Item  3 

Reference  and  geometric  information. 

LROUT  U  f) 

Reference  length  used  in  moment  calculations. 

SROUT  (Sref) 

Reference  area  used  in  force  and  moment  calculations 

XMC  (x  ) 
m 

Moment  center  of  missile  measured  from  missile  sta¬ 
tion  zero,  dimensional. 

A  (a) 

Radius  of  missile,  dimensional. 

TIPRAD 

Nose  tip  radius,  dimensional. 
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ETAN  (n) 

Half  angle  of  body  nose  for  pointed  body;  or  angle 
between  tangent  to  nose  at  juncture  of  spherical 
cap  and  rest  of  nose  and  body  axis,  degrees. 

XCLE 

Distance  from  missile  station  zero  to  leading  edge 
of  canard  root  chord,  dimensional. 

Item  4 


CMACH(J) , 
(Moo(J)) 

Free-stream  Mach  number;  1  <  J  <  NMACH . 

Item  5 


XNOSE(M) 

Axial  location  entries  in  nose  coordinate  table; 

1  <.  M  <.  NNOSE ;  XNOSE(l)  is  axial  location  of  nose 
tip  from  MS  0,  dimensional. 

RNOSE(M) 

Corresponding  radial  location  entries  in  nose  coor¬ 
dinate  table;  1  <  M  <  NNOSE,  dimensional. 

Item  7 


ALFAC(K) , 
(ac(K)  ) 

Body  angle  of  attack  in  degrees;  1  <  K  <  NALFA. 

Item  8 


PHI(L) , 

( 4>  (L)  ) 

Bank  angle  in  degrees;  angle  between  z  and  zQ  axes; 
positive  measured  clockwise  viewed  from  rear; 

1  <  L  <  NPHI . 

Item  9 


DEFC ( I , J) , 

<6C(I)) 

Jth  set  of  canard  fin  deflection  angles,  degrees. 

1  <  I  <  4;  number  of  the  canard  fin. 

1  <  J  <  NDELT . 

Item  10 

Canard  geometry. 

SPANC ( s  ) 
m 

Maximum  semispan,  measured  from  body  centerline, 
of  canard  fins,  dimensional. 

XCHL  (xRL) 

Axial  distance  to  canard  hinge  line,  meausred  from 
MS  0,  dimensional. 

XCTIPL 

Axial  distance  to  leading  edge  of  canard  tip  chord, 
measured  from  MS  0,  dimensional;  XCTIPL  >  XCLE. 

XCTE 

Axial  distance  to  canard  trailing  edge,  measured 
from  MS  0,  dimensional. 
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Item  11 


Fin  airfoil  data 


NFOILS 


Number  of  strips  across  span  which  the  fin  is  to 
be  divided  into;  1  <  NFOILS  <  10. 


TOC  (t/c)  Ratio  of  flat-plate  thickness  to  airfoil  chord 

length . 


FOC  (f/c)  Ratio  of  flat-plate  length  to  airfoil  chord  length. 


Item  12 


999 


This  card  causes  the  program  to  stop  searching  for 
more  data  and  the  run  is  stopped. 


A. 6  SAMPLE  CASE  INPUT  DATA 


The  input  data  for  a  sample  case  is  presented  in  this  section 
to  illustrate  the  use  of  the  computer  program.  The  configuration 
used  is  shown  in  figures  A. 5  and  A. 6.  The  body  used  is  NIB  from 
Reference  A7  but  is  terminated  at  the  canard  trailing  edge.  The 
canard  fin  is  one  of  the  fins  used  in  the  parametric  study  for  the 
present  report. 

The  input  data  deck  for  this  case  is  shown  in  figure  A. 7. 

The  calculation  is  done  for  five  Mach  numbers,  five  sets  of  fin 
deflection  angles,  five  angles  of  attack,  and  one  roll  angle. 
Eleven  points  are  used  in  the  Simpson's  rule  integrations  in  sub¬ 
routine  REVFLO  and  there  are  16  entries  in  the  nose  shape  table. 
The  nose  normal-force-curve  slope  is  not  input  (NCA  =  0)  and  the 
influence  of  nose  vortices  is  calculated  to  the  canard  root-chord 
leading  edge  (NVORT  =  1) .  The  option  of  printing  only  the  fin 
loads  and  the  maximum  and  minimum  hinge-moment  tables  is  selected 
(NOUTPT  =  0) .  The  crossflow  is  considered  turbulent  (TURB  =  T) 
and  the  input  reference  area  and  reference  length  are  used  for  the 
fin  loads  (REFER  =  F) . 


A. 7  DESCRIPTION  OF  OUTPUT 


This  section  of  the  report  will  describe  the  output  from  the 
program.  The  output  for  the  sample  case  will  be  used  for  this 
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purpose.  This  is  shown  in  figure  A. 8.  Not  all  of  the  output  is 
included  in  the  figure.  Sufficient  is  presented  to  describe 
everything  and  also  to  verify  that  the  program  is  working  properly. 

The  first  section  of  output  is  the  input  data  and  is  shown 
in  figure  A. 8 (a).  The  user  provided  heading  information  is  first 
printed.  This  is  followed  by  the  missile  geometry  and  the  refer¬ 
ence  lengths  and  areas.  The  quantities  CPXRC  and  CPYRC  are  the 
canard  root  chord  and  exposed  span,  respectively.  The  remaining 
output  in  the  figure  are  the  flight  condition  parameters  and 
program  operation  indices.  All  of  this  is  printed  for  both  output 
options,  NOUTPT  =  0  or  1,  at  the  beginning  of  each  case. 

The  second  section  of  output  is  shown  in  figure  A. 8(b).  Here 
the  fin-alone  normal-force  curve  and  center-of-pressure  curve  are 
tabulated  for  a  particular  Mach  number.  This  information  is 
printed  each  time  the  Mach  number  changes.  The  normal-force  curve 
is  obtained  from  the  data  base  in  the  program.  The  center-of- 
pressure  location  as  a  fraction  of  the  root  chord  (CR)  is  calcu¬ 
lated  in  subroutine  XCPVSA.  For  this  case  the  leading-edge  shock 
wave  is  detached  so  the  Busemann  theory  is  used  for  the  thickness 
correction  at  a  =  2°.  If  the  shock  is  attached,  shock-expansion 
theory  is  used  at  the  low  angles  of  attack  and  is  so  indicated  in 
the  output.  The  a  =  45°  value  is  obtained  from  the  data  base  in 
the  program.  This  output  is  printed  for  both  output  options. 

Figure  A. 8(c)  shows  the  canard  fin  force  and  moment  output 
for  NOUTPT  =0.  At  the  top  of  the  page  the  Mach  number,  canard 
deflection  angles,  and  roll  angle  are  printed.  This  is  followed 
by  the  canard  fin  forces  and  moments  for  each  of  the  angles  of 
attack,  in  this  case  five  angles.  For  each  angle  of  attack,  the 
canard-fin  nose-vortex-induced  equivalent  angle  of  attack  (DALFEQ)V 
is  tabulated  as  are  the  vortex  induced  normal-force  coefficient 
(DCN)V,  root-bending-moment  coefficient  (DCBM)V,  hinge-moment 
coefficient  (DCHM)V,  and  rolling-moment  coefficient  (DCRM)V. 

These  are  followed  by  the  total  canard-fin  equivalent  angle  of 
attack  ALFEQ  and  the  normal-force  CNF,  root-bending-moment  CBM, 
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hinge-moment  CHM,  and  rolling-moment  CRM  coefficients.  The  last 
column  lists  the  fin  center-of-pressure  location  as  a  fraction 
of  the  root  chord.  All  of  the  coefficients  use  SROUT  and  LROUT, 
input  quantities,  as  the  reference  area  and  reference  length. 

The  positive  directions  of  the  force  and  moment  coefficients  are 
shown  in  figure  A. 9.  A  page  of  output  like  this  is  printed  for 
each  Mach  number,  fin-deflection  angle,  and  roll  angle  combination 
for  output  option  NOUTPT  =  0,  the  one  used  in  the  sample  case. 

When  the  NOUTPT  =  1  output  option  is  used,  the  output  for  each 
angle  of  attack  is  expanded.  An  example  for  the  a  =  20°  point 
shown  in  figure  A. 8(c)  is  presented  in  figure  A. 8(d) .  Output  like 
this  is  printed  for  each  Mach  number,  canard  deflection  angle, 
roll  angle,  and  angle  of  attack  combination.  For  the  sample  case 
this  would  be  125  pages  of  output  like  that  shown  in  figure  A. 8(d) . 

At  the  top  of  the  page  the  Mach  number,  canard  deflection 
angles,  and  roll  angle  are  listed.  This  is  followed  by  the  re¬ 
sults  for  the  nose  section.  The  strengths  and  positions  of  the 
nose  vortices  at  the  canard  root-chord  leading  edge  are  tabulated. 
The  strengths  are  nondimensionalized  by  2i\Voa  and  the  positions 
in  the  xq,  y  ,  zq  coordinate  system  by  the  body  radius,  a.  The 
coordinate  system  is  shown  in  figure  A. 10.  The  last  nose  sec¬ 
tion  data  tabulated  are  the  forces  and  moments  in  the  unrolled 
(x  ,  y  ,  z  )  coordinate  system.  Positive  directions  of  the 
forces  and  moments  are  shown  in  figure  A. 10. 

The  first  output  printed  for  the  canard  section  is  that  which 
was  tabulated  for  the  NOUTPT  =  0  option,  the  fin  loads  made  dimen¬ 
sionless  by  SROUT  and  LROUT.  This  is  followed  by  a  listing  of 
fin  hinge-moment  coefficient  as  a  function  of  hinge-line  location. 
The  hinge-line  location  is  listed  as  a  fraction  of  the  fin  root- 
chord  length.  The  next  output  again  lists  the  fin  forces  and 
moments  but  this  time  the  nondimensionalizing  quantities  are  the 
fin  planform  area,  root  chord,  and  exposed  semispan.  The  final 
output  tabulated  for  the  canard  section  are  the  canard  section 
contributions  to  the  total  loads.  These  include  body  loads. 
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The  forces  and  moments  are  tabulated  in  the  unrolled  coordinate 
system  (xQ,  yQ,  zQ)  and  rolled,  or  body-axis,  system  (x,  y,  z) 
shown  in  figure  A. 10. 

The  final  output  tabulated  in  figure  A. 8(d)  is  a  summary  of 
the  total  loads  on  the  configuration.  The  nose  loads,  canard 
section  loads,  and  the  sum  of  these  two  are  listed. 

The  last  output  printed  for  both  output  options  is  shown  in 
figure  A. 8(e).  This  is  a  tabulation  of  the  maximum  and  minimum 
hinge-moment  coefficients  for  fin  4  as  a  function  of  hinge-line 
position  and  Mach  number.  A  table  like  this  is  printed  for  each 
roll  angle.  The  values  listed  are  obtained,  for  each  hinge-line 
location  and  Mach  number,  by  searching  through  all  of  the  angle 
of  attack  and  fin-deflection  angle  combinations  for  which  calcu¬ 
lations  were  made. 
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Figure  A. 2.-  Flow  of  program  computations. 
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Common  blocks  in  program. 
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Figure  A. 4.-  Continued. 


Item  8  Format (BF10. 5) ,  8  values  of  PHI  per  card 
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Figure  A, 5-  Configuration  used  for  sample  case. 
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Figure  A. 6.-  Dimensions  of  canard  fins. 
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Figure  A. 7.-  Input  data  for  sample  case. 
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2 

3 


(DALFEQ) V 
0.00 
0.00 
0.00 
0.00 


1DCN1 V 
0.0000 
0.0000 
0.0000 
0.0000 


(DCBM)V 

0.0000 

0.0000 

0.0000 

0.0000 


IDCHM) v 
0.0000 
0.0000 
0.0000 
0.0000 


(DCRM) V 
0.0000 
0.0000 
0.0000 
0.0000 


FIN  LOADS  * 

ON 

SROUT  ANO  LROUT) 

CRM 

XCPCR 

ALFEQ  CNF 

CBM 

CHM 

• 

-.00  -.0000 

-.0000 

-.0000 

-.0100 

•  6136 

• 

22,75  .7620 

.2492 

.0011 

•  6402 

.6310 

• 

-.00  -.0000 

-.0000 

-.0000 

-.0000 

.6136 

• 

22.75  .7820 

.2492 

.0011 

.6402 

.6310 

effect 

OF  HINGE-LINE  LOCATION  ON  HINGE 
HINGE  MOMENTS 

MOMENT 

1L/CR 

FIN  1 

FIN  2 

FIN  3 

FIN  4 

.2000 

•  0000 

-.4128 

.  0000 

-.4126 

.3000 

.0000 

-.3170 

.0000 

-.3170 

.4000 

•  0000 

-.2213 

.0000 

-.2213 

.5000 

.0000 

-.1255 

.0000 

-.1255 

.6000 

.0000 

-.0297 

.0000 

-.0297 

.7000 

-.0000 

.0661 

-.0000 

.0661 

.6000 

-.0000 

.1616 

-.0000 

.1618 

total  CANARO  fin  LOADS  based  ON  PLANFORM  area,  root  choro.  ano  exposed  SEMISPAN 


panel 

1 

2 

3 


ALFEQ 

-.00 

22.75 

-.00 

22.75 


CNF 

-.0000 

1.0919 

-.0000 

1.0919 


unrolled  COORDINATES 


CAO  ■ 
CYO  • 

czo  « 

CL  ■ 
CD  - 


O.OOOrt 
.0000 
1.950? 
1.8325 
•  6670 


Cm  aO 
CMYO 
CM20 


0.0000 

-11.6539 

.0000 


CBM 

CHM 

CRM 

-.0000 

-.0000 

-.0000 

.3768 

.0013 

.9731 

-.0000 

-.0000 

-.0000 

.3788 

.0013 

.9731 

TO  TOTAL  LOADS 

•  •• 

AND  LROUT) 

ROLLED 

COORDINATES 

cx  ■ 

0.0000 

C  M  X  B  0.0000 

CY  » 

.0000 

CMY  B  -11*6539 

CZ  ■ 

1 .9502 

CM?  b  .0000 

•••«  sumnahy  Of 

ALPHA  B  20.00 


total  loads  ••• 

PHI  ■  0.00 


UNROLLFO  COORDINATES 


NOSF 

CANARD 

totals 


cxo 

0.00000 

0.00000 

0.00000 


CYO 

.00000 

.oooon 

.00000 


C70 
1 .30596 
1.95015 
3.25611 


CMXO 
0.00000 
0.00000 
0. oooon 


CMYO 
-2.B4012 
-11.65393 
-14.49405 


•  •• 

CMZO  c* 

0.00000  o.ooooo 
,00000  0.00000 
.00000  0.00000 


BODY  AXIS  COORDINATES 


CY 

0.00000 

.00000 

.00000 


c? 

1.30596 

1.95015 

3.25611 


CM  X  CMY 

0.00000  -2.6401? 

0.00000  -11.65393 
0.00000  -14.49a  05 


CM? 

0.00000 

•00000 

.00000 


A X  ]  A!  CENTERS  Of  PRESSURE 


NORMAL  FORCE  (CPX) 
SIDE  FO»CE  (CPY) 


4.45134 

5.9604? 


(d)  Forces  and  moments  for  NOUTPT  =  1 
Figure  A. 8.-  (Continued). 
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MAXIMUM  AND  minimum  HINGE  moments  FQR  fin  4  (BASED  ON  SROUT  AND  LROUT ) .  OVER  THE  RANGES  OF  ANGLE  OF  ATTACK 
AND  FIN  DEFLFCTION  ANGLE  FOR  WHICH  CALCULATIONS  WERE  MAUL,  I- OR  A  GIVEN  ROLL  ANGLE.  THEY  ARE  TABULATED 

as  a  function  of  hinge-line  location  and  mach  number. 

ROLL  ANGLE  (PHI)  =  0.00 

•••*•••••*••«  MAXIMUM  hinge  MOMENT  ***** . . . MINIMUM  HINGE  MOMENT . ******** 


xhl/cr 

M* 1.20 

M=1 .60 

.?o 

.3461 

.3801 

.21 

.3379 

.3713 

.22 

.3298 

.3626 

.23 

.3216 

.3536 

.24 

.3135 

.3450 

.25 

.3054 

.3363 

.26 

.2972 

.3275 

.27 

.2891 

.3180 

.28 

.2809 

.3100 

.29 

.2728 

.3012 

.30 

.2647 

.2925 

•  31 

.2565 

.2837 

.32 

.2484 

.2749 

.33 

.2402 

.2662 

.34 

.2321 

.2574 

.35 

.2239 

.2486 

.36 

.2158 

.2399 

.37 

.2077 

.2311 

.38 

.1995 

.2224 

.39 

.1914 

.2136 

.40 

.1832 

.2048 

.41 

.1751 

.1961 

.42 

.1670 

.1873 

.43 

.1588 

.1785 

•  44 

.1507 

.1690 

.45 

.1425 

.1610 

.46 

.1344 

.1522 

.47 

.1263 

.1435 

•  48 

.1181 

.1347 

.49 

.1100 

.1260 

.50 

.1016 

.1172 

.51 

.0937 

.1084 

.52 

.0856 

.0997 

.53 

.0774 

.0909 

.54 

.0693 

.0821 

.55 

.0611 

.0734 

.56 

.0530 

.0646 

.57 

.0449 

.0558 

.58 

.0367 

.0471 

.59 

.0286 

.0383 

.60 

.0204 

.0296 

•  61 

.0123 

.0208 

.62 

.004? 

.0120 

.63 

.0045 

.0061 

•  64 

.0108 

.0208 

.65 

.0186 

.0354 

.66 

.0278 

.0501 

.67 

.0377 

.  .0648 

■  68 

.0483 

.0794 

.69 

.0596 

.094  1 

.70 

.0716 

.  1  088 

.71 

•  0841 

.1234 

.72 

.0969 

•  13«1 

.73 

.  1097 

.1528 

.7* 

.  1226 

.  1674 

.75 

.1354 

.1821 

.76 

.148? 

.  1  9*8 

.  77 

.1611 

.2114 

•  78 

.1739 

.2261 

.79 

.1867 

.2408 

.80 

.  1996 

.2554 

M=2 .00 

M  =  2  »  50 

M=  3 .00 

.4078 

.4529 

•  4944 

.3984 

.4424 

.4029 

.3890 

.4320 

.4714 

.3796 

.4216 

.4599 

.370? 

.4111 

.4484 

.3608 

.4007 

.4369 

.3515 

.3902 

.4254 

.3421 

.3790 

.4139 

.3327 

.3693 

.4024 

.3233 

.  3589 

.3909 

.3139 

.3484 

.3793 

.3045 

.3380 

.3678 

.2951 

.3275 

.3563 

.2857 

.3171 

.3448 

.2764 

.3066 

.3333 

.2670 

.2962 

.3218 

.2576 

.2857 

.3103 

.2482 

.2753 

.2988 

.2388 

.2648 

.2873 

.2294 

.2544 

.2758 

.2200 

.2440 

.2643 

.2106 

.2335 

.2528 

.2013 

.2231 

.2413 

.1919 

.2126 

.2298 

.1825 

.2022 

.2182 

.1731 

.1917 

.2067 

.  1637 

.1813 

.1952 

.1543 

.  1700 

.1837 

.1449 

.  1604 

.1722 

.1356 

.  1499 

.1607 

.  1262 

.  1  395 

.1492 

.  1  168 

.1290 

.1377 

.1074 

.1186 

.  1262 

.0980 

.1081 

.  1147 

.0886 

.09/7 

.1032 

.079? 

.0873 

.0917 

.0698 

.0768 

.080? 

.  0605 

.0664 

.0687 

.0511 

.0559 

.0571 

.0417 

.0455 

.0456 

.0323 

.0350 

.0341 

.0229 

.0246 

.0254 

.0140 

.0158 

.017? 

.0113 

.0085 

.0089 

.0280 

.0226 

.0090 

.0447 

.0419 

.0306 

.0614 

.0612 

.0523 

.  0781 

.0804 

.0739 

.0949 

.0997 

.0956 

.1116 

.1190 

.1172 

.1283 

.1302 

.1309 

.1450 

.  1575 

.1606 

.1617 

.  1  768 

•  1821 

.  1  785 

.  1960 

.2038 

.  195? 

.2153 

.2264 

.2119 

.2346 

.2471 

.2286 

.2538 

.268  1 

.2453 

.  ? 73  1 

.  2904 

.262] 

.2924 

.3120 

.2788 

.3117 

.3337 

.2955 

.  3309 

.3553 

XHL/CR 

M=1 .20 

M= 1 .60 

.20 

-.5706 

-.6246 

.21 

-.6577 

-.6099 

.22 

-.5449 

-.5953 

.23 

-.5321 

-.5806 

.24 

-.519? 

-.5659 

.25 

-.5064 

-.5513 

.26 

-.4936 

-.5366 

.27 

-.4807 

-.5219 

.28 

-.4679 

-.5073 

.29 

-.4551 

-.4926 

.30 

-.4422 

-.4779 

.31 

-.4294 

-.4633 

.32 

-.4165 

-.4486 

.33 

-.4037 

-.4339 

.34 

-.3909 

-.4193 

.35 

-.3780 

-.4046 

.36 

-.3652 

-.3899 

.37 

-.3524 

-.3753 

.38 

-.3395 

-.3606 

.39 

-  •  3?67 

-.3459 

.40 

-.3139 

-.3313 

.41 

-.3010 

-.3166 

.42 

-.2882 

-.3019 

.43 

-.2754 

-.2873 

.44 

-.2625 

-.2726 

.45 

-.2497 

-.2579 

.46 

-.2368 

-.2433 

.4  7 

-.2240 

-.2286 

.48 

-.2112 

-.2139 

.49 

-.1983 

-.1993 

.50 

-.1855 

-.  1846 

.51 

-.1727 

-.1699 

.52 

-. 1598 

-.1553 

.53 

-.1470 

-.1406 

.54 

-.1342 

-.1259 

.55 

-.1213 

-.1113 

.56 

-.1085 

-.0966 

.57 

-.0956 

-.0819 

.58 

-.0828 

-.0673 

.59 

-.0700 

-.0526 

.60 

-.0571 

-.0383 

.61 

-.0443 

-.0248 

.62 

-.0315 

-.0129 

.63 

-.0186 

-.0039 

.64 

-.0121 

-.0055 

.65 

- • 0  ?03 

-.0U3 

.66 

-.0284 

-.0230 

.67 

-.0365 

-.0318 

.68 

-.0447 

-.0406 

.69 

-.0528 

-.0493 

.70 

-.0610 

-.0501 

.71 

-.0691 

-.0668 

.72 

-.0772 

-.0756 

.73 

-.0854 

-.0844 

.74 

-.0935 

-.0931 

.76 

-.1017 

-.1019 

.76 

-.1 090 

-.1107 

.77 

-.1180 

-.1194 

.  78 

-. 1?61 

-.1282 

.79 

-.1 34? 

-.1370 

.80 

-.1424 

-. 1457 

M  =  2*  00 

M=? • 50 

M  =  3 • 00 

-.7077 

-.8252 

-.9436 

-.6910 

-.0059 

-.9219 

-.6743 

-.7867 

-.9003 

-.6575 

-.7674 

-.8786 

-.6408 

-.7481 

-.8570 

-.6241 

-.7289 

-.8353 

-.6074 

-. 7096 

-.8137 

-.5907 

-.6903 

-.7920 

-.5739 

-.6711 

-.7704 

-.5572 

-.6518 

-.7487 

-.5405 

-.6325 

-.7271 

-.5238 

-.6133 

-.7054 

-.5071 

-.5940 

-.6838 

-.4903 

-.5747 

-.6621 

-.4736 

-.5554 

-.6405 

-.4569 

-.5362 

-.6188 

-.4402 

-.5169 

-.5972 

-.4235 

-.4976 

-.5755 

-.4067 

-.4704 

-.5539 

-.3900 

-.4591 

-.5322 

-.3733 

-.4398 

-.5106 

-.3566 

-.4206 

-.4809 

-.3399 

-.4013 

-.4673 

-.3231 

-.3820 

-.4457 

-.3064 

-.3628 

-.4240 

-.2897 

-.3435 

-.4024 

-.2730 

-.3242 

-.3807 

-.2563 

-.3050 

-.3591 

-.2395 

-.2857 

-.3374 

-.2228 

-.2664 

-.3158 

-.2061 

-.2471 

-.2941 

-.1894 

-.2279 

-.2725 

-.1727 

-.2086 

-.2508 

-.1559 

-.1893 

-.2292 

-.139? 

-.1701 

-.2075 

-.1225 

-.1508 

-.1659 

-.1058 

-.1315 

-.1642 

-.0891 

-.1123 

-.  1426 

-.0734 

-.0943 

-.1209 

-.0577 

-.0763 

-.0995 

-.0434 

-.0596 

-.0795 

-.0300 

-.0439 

-.0595 

-.0182 

-.0295 

-.0394 

-.0091 

-.0156 

-.0194 

-.0052 

-.0068 

-.0119 

-.0146 

-.0172 

-.0234 

-.0240 

-.0277 

-.0349 

-.0334 

-.0381 

-.0464 

-.0428 

-.0486 

-.0579 

-.0522 

-.0590 

-.0694 

-.0616 

-.0694 

-.0809 

-.0710 

-.0799 

-.0925 

-.0003 

-.0903 

-.1040 

-.0897 

-.  1  008 

-.1155 

-.0991 

-.1112 

-.1270 

-.1005 

-.1217 

-.1385 

-.1179 

-.1 321 

-.1500 

-.1273 

1426 

1615 

-.1367 

-.1530 

-.1730 

-.1461 

-.1635 

-.1845 

-.1554 

-. 1739 

-.  1960 

END  OF  CALCULATIONS  FOR  THIS  CASE 

(e)  Summary  of  maximum  and  mimimum  hinge  moments 
Figure  A. 8.-  (Concluded). 
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CBM1,  CRM1 


CBM2 

CRM2 


Hinge 

line 


CBMi 

CRMi 


i 


Figure  A. 9.-  Positive  directions  for  canard  fin 
force  and  moment  coefficients. 
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Figure  A. 10- 
coordinate 


Unrolled  (xq,  y  ,  z  ,)  and  rolled  (x,  y,  z) 
systems  and  positive  directions  of  forces 
and  moments  in  the  two  systems.. 
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